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ABSTRACT 


The  physical  properties  and  thermochemistry  of  decomposition  of  Si05  at 
high  temperatures  has  been  thoroughly  reviewed  and  detailed  calculations 
have  been  made  of  thermodynamic  properties  of  interest  for  re-entry  vehicle 
applications.  ./Available  thermodynamic  data  indicate  that  the  primary  reac¬ 
tion  involved  in  the  vaporization  of  si02  under  neutral  conditions  is 


Sio2(l)  — sio(g)  ♦  ^02(f) 

The  vapor  pressure  relation  corresponding  to  this  equilibrium  is 
Pt  -  3/2 


(1) 


0--^) 

e  1 


(2) 


T  iii°talJ?P°r  pre,,ure  in  atmosphere.  A  normal  boiling  point  of 
3° 8 5  K  is  predicted  from  this  relation.  If  other  equilibria  are  considered, 
including  formation  of  SiO^^ .  dimer  formation,  and  dissociation  of  oxygen. 

a  more  refined  boiling  point  of  3070*  ±  70*  K  (28 00*  C)  is  obtained.  This 
value  has  been  compared  with  data  in  the  literature  and  is  consistent  with  the 
most  reliable  data  available.  When  Si02  i.  heated  from  room  temperature 
to  its  boiling  point  and  decomposed  at  that  temperature,  the  total  heat  absorp¬ 
tion  is  calculated  to  be  6650  Btu/lb.  F 


Many  other  physical  properties  of  si02  which  are  of  interest  for  re-entry 

vehicle  design  purposes,  have  been  reviewed  and  are  summarized  in  an 
appendix  at  the  end  of  the  report. 
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I.  INTRODUCTION 


One  of  the  most  promising  inorganic  compound”  for  potential  use  as  a  con¬ 
sumable  material  on  ablating  nose  cones  is  fused  silica.  *  To  predict  the 
behavior  of  this  material  when  subjected  to  the  high  heating  rates  of  re-entry, 
it  is  necessary  to  know  much  about  its  thermochemistry  of  decomposition  at 
high  temperatures.  Once  the  vapor  pressure  of  silica  under  equilibrium 
conditions  is  known,  then  it  may  be  possible  to  predict  its  vaporization  be¬ 
havior  under  the  nonequilibrium  conditions  of  re-entry.  This  report  presents 
an  extensive  analysis  of  the  thermochemistry  of  silica.  Other  properties 
such  as  viscosity  of  the  melt,  thermal  shock  resistance,  etc.  are  also  of 
interest  and  have  been  included. 

Although  the  glass  industry  is  an  old  one  and  SiC>2  is  one  of  the  most  plenti¬ 
ful  compounds  on  the  earth's  crust,  there  exists  a  paucity  of  accurate  vapor 
pressure  data  for  pure  si02  at  elevated  temperatures.  The  data  of  Brewer 
and  Mastick,^  and  Porter,  Chupka,  and  Inghram  ^  represent  the  most  recent 
available  data  on  the  vapor  pressure  of  S1O2  .  Their  data  were  all  obtained 
below  2,  000"  K  but  in  this  memorandum  they  are  shown  to  agree  satisfactorily 
with  vapor  pressures  calculated  from  other  available  thermodynamic  data. 

The  tabulated  boiling  points  for  siCK  ,  which  are  recorded  in  various  hand¬ 
books,  have  values  of  2,  230*  C,  2,  590°  C  and  2,  950*  C.  These  have  all  been 
traced  to  the  work  of  Ruff  and  co-authors.  All  of  boiling  point 

determinations  can  be  criticized  because  they  were  carried  out  in  a  reducing 
atmosphere.  Since  a  reducing  atmosphere  leads  to  a  higher  vapor  pressure 
than  would  be  expected,  "the  lower  boiling  points  are  not  considered  to  represent 
the  properties  of  the  material.  Moreover,  the  two  higher  boiling  points, 

2,  590°  C  and  2,950°  C,  have  been  based  on  a  single  vapor  pressure  determina¬ 
tion  --  11  mm  at  2,060°  C.  In  their  original  paper,  Ruff  et  al  ^  stated  that 
this  measurement  yields  2,  590°  C  as  a  boiling  point.  Some  nine  years  later 
Ruff  stated  without  further  evidence  that  this  was  in  error  and  should  be 
2,  950°  C.  A  review  of  Ruff's  extrapolation  has  indicated  that  he  was  probably 
correct  the  first  time  in  evaluating  his  dat'  to  obtain  a  normal  boiling  point 
of  2,  590°  C.  However,  the  original  criticism  that  his  measurement  was  done 
under  reducing  conditions  is  still  valid,  and  the  actual  boiling  point  is  prob¬ 
ably  above  2,  590°  C  as  he  himself  stated  in  his  original  paper.  ^ 

In  the  present  paper  various  possible  reactions,  which  can  occur  when  SiC^ 
vaporizes,  have  been  considered.  From  available  thermochemical  data  the 
free  energy  change  for  each  reaction  has  been  obtained  as  a  function  of 
temperature.  Hence,  equilibrium  constants  for  these  reactions  can  also  be 
given  as  functions  of  temperature.  It  is  then  possible  to  estimate  the  various 
species  or  products  present  at  equilibrium  at  different  temperatures  and 
pressures. 
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II.  THERMOCHEMISTRY  OF  THE  REACTION 


The  above  reaction  has  been  considered  by  Brewer  et  al8*  2  and  Porter  et  al  3 
to  be  the  most  important  one  occurring  when  Si02  decomposes  under  neutral 
conditions.  The  results  of  the  present  study  are  consistent  with  these  earlier 
findings.  Under  oxidizing  conditions  (i.  e.  ,  with  excess  oxygen  present),  the 
above  reaction  can  be  shifted  to  the  left  so  that  S1O2  gas  formation  is  re¬ 
pressed. 


The  equilibrium  constant  of  the  reaction  and  the  equilibrium  vapor  pressure 
of  liquid  SI02  can  be  obtained  by  calculation  from  the  heat  of  reaction  and 
entropy  change  of  the  reaction.  First  we  calculate  the  heat  of  reaction  and 
entropy  change  for  the  reaction  from  available  data. 

A.  HEAT  OF  REACTION 


To  calculate  the  heat  of  the  above  reaction  we  can  subtract  the  heat  of  forma¬ 
tion  of  S1O20)  from  the  heat  of  formation  of  the  products.  Since  the  heat  of 

formation  of  the  elements  is  taken  to  be  zero,  it  is  necessary  to  evaluate  the 
heat  of  formation  of  Si02^j  and  SiO(g)  only.  For  the  present  case,  the 

temperature  range  from  2,  000  to  3,  000*  K  has  been  considered  of  primary 
interest.  Because  data  for  this  temperature  range  are  not  available,  it  is 
necessary  to  use  existing  data  for  lower  temperatures  and  calculate  the 
thermodynamic  functions  at  the  higher  temperatures  by  extrapolation. 

1.  Heat  of  Formation  of  SIO-, 

- 2(1) 

9 

Coughlin  has  tabulated  values  for  the  heat  and  free  energy  of  formation 
of  Si02  up  to  2,  000*  K.  His  results  are  based  on  the  heat  of  formation 
data  of  Humphrey  and  King.  10  His  heat  of  formation  data  at  2,  000*  K  are 


^(l)  ♦  %>^Si%)  ^2,000*  "  "215,200  ±  1,000 


SiO)  *  °2(„— •  SW2(0c,i.«.b.li«e)AH*.OOO<>C  *  ±  250  (4) 

/9-crietobalite  has  a  melting  point  of  1,723*  C  (1,996*  K).  according  to 
Levin  et  al,  and  a  value  of  1, 728*  C  (2,  001*  K),  according  to  Campbell11 
For  convenience,  we  can  take  2,000*K  as  the  melting  point  of 
cristobalite  with  negligible  error.  The  diagrams  of  Levin  et  al12 


indicate  that  silica  glass  is  the  stable  condensed  species  of  Si<>2 
above  2,  000*  K.  In  further  discussions,  the  terms  silica  glass  or  silica 
liquid  are  used  interchangeably  for  the  range  above  2,  000*  K. 

la  any  thermochemical  calculation,  it  is  desirable  to  start  with  the  most 
accurate  data  available.  Referring  to  the  reactions  (3)  and  (4),  it  might 
seem  preferable  to  start  with  reaction  (4)  rather  than  reaction  (3)  since 
its  heat  of  formation  is  given  to  a  greater  precision.  However,  the  un¬ 
certainty  in  the  heat  of  fusion  of  /9  -cristobalite  makes  either  reaction  a 
suitable  starting  point.  Thus,  we  take  the  difference  between  reactions 
(4)  and  (3)  as  the  heat  of  fusion  of  0  -cristobalite,  AHf  •  2500  calories/mole  • 
From  the  entropy  data  of  Kelley13*  14we  calculate  a  heat  of  fusion  at 
2,  000*  K  of  2,  060  calories /mole.  A  value  of  3,  600  calories /mole  is 
quoted  by  Kubaschewski  and  Evans15.  The  NBS  compilation 16  gives 
a  value  of  1,840  calories /mole.  The  original  data  quoted  by  Kelley1? 
indicates  experimental  values  for  the  heat  of  fusion  of  cristobalite 
to  be  in  the  range  of  1,000  to  3,  700  calories /mole.  He  derives 
2,  100  calories /mole  as  an  average  value.  This  range  in  values  for  the 
heat  of  fusion  shows  why  the  data  for  reaction  (3)  are  less  precise,  but 
•till  the  best  available.  Therefore,  the  data  for  reaction  (3)  are  used 
for  the  present  calculations. 

The  heat  of  formation  at  temperatures  above  2,  000*  K  can  be  calculated 
by  means  of  the  Kirchoff  relation: 

T 

*  ^2000  ♦  f  ACp  dT 
Or  *2000 


-  ^2000  +  <HT  -  HWsiOj^  ~<Ht  -H2000>sia)  —  (Rj.  -  (5) 

The  results  of  such  a  calculation  at  2,  000*  .  2,  500*  ,  and  3,  000*  K  are 
given  in  Table  I.  For  this  calculation,  the  heat  content  of  SiOjd)  was 
estimated  by  assuming  that  its  average  heat  capacity  was  1  calorie  /mole 
degree  greater  than  that  given  by  Kelley13  for  at  2,000*  K. 

That  is,  C  -  21.66  was  assumed. 

P 

According  to  Kubaschewski  and  Evans15,  a  rough  estimate  of  the 
heat  capacity  of  molten  liquids  can  be  obtained  by  assuming  that 
each  atom  contributes  7.  25  calories/g.  atom*K.  This  leads  to  a  value 
°£  21.  75  calorie/mole  *  K.  In  the  paper  by  Humphrey  et^al1®,  the  heat 
capacity  of  liquid  silica  was  taken  to  be  Cp  —  24.0  •  In  the  present  case. 
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4  ®  extraP°latlon  to  High-temperatures  leads  to  considerable  uncertainty 
in  the  heat  content.  The  uncertainty  in  the  heat  capacity  is  roughly 
estimated  to  be  1  calorie  /mol*  degree  thus  giving  CpsiO^  -  21.66  ±1.0  . 

Sources  of  data  for  the  calculations  of  Table  I  are  given  in  the  table. 
Rough  estimates  of  the  uncertainties  in  the  data  have  been  made.  The 

f°r  thC  heat  °f  formation  the  temperature  range  2,  000  to 
.5.000  K  become  more  uncertain  as  temperature  rises  because  of  the 


uncertainty  in  the  heat  capacity  of  liquid  si03  in  thi 


s  range. 


2. 


Heat  of  Formation  of  SiO^j 


The  heat  of  formation  of  SiO(g) 
calories.  This  value  is  the  one 
It  is  based  on  a  recalculation  of 
studied  the  reaction 


at  298*  K  is  taken  to  be  -2 1 , 41 1  ±  574 
recommended  by  Edwards  and  Brewer.®*  ^ 
the  work  of  Schafer  and  Hornle20who  had 


*<•>  *  2Si°(s)  * 

below’ *°urcc*  of  con»idered  by  Brewer  and  Edwards,  8  are  given 

Investigators  and  Equilibrium  Studied  A»P298  (K  calorie  /mole) 


Grube  and  Speidel21 

"*«>  *  H2%)  +  *°(l> 

Brewer  and  Mastick2 

-15,  291  ±  2,  000 

*%>—  *<«>  +  *%> 

Gel'd  and  Kochnev 

-21,  159  *  2,  000 

*<•>  ♦  S*0*.)—*  2Si°(g) 

Schafer  and  Hornle  20 

-18,  212  ±  960 

*<«)  ♦  SiO^^  — *  2SiO(f) 

Tombs  and  Welch23 

-21,411  ±  570 

"*,)  +  Si°2(.)—  H2°(«)  ♦  “Of,) 

Low  temperature 
-26,  004  ±  2,  400 

High  temperature  data 
-21, 269  ±  1.  500 
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Brewer  end  Edwards  considered  each  of  the  above  sources  very  carefully 
and  chose  the  data  of  Schafer  and  Hornle  as  being  the  best.  It  should  be 
pointed  out  that  the  arguments  they  use,  based  on  dissociation  energies 
of  SiO^  ,  were  fortuitously  correct  inasmuch  as  more  recent  data  shows 
they  used  too  low  a  value  for  the  heat  of  sublimation  of  silicon  (88.  9 
Kcalorie/mole)  and  for  the  dissociation  energy  of  sio(g)  ,  (169  Kcalorie/ 
mole).  A  more  recent  value  for  the  heat  of  sublimation  of  monatomic 
silicon  is  105  Kcalorie/mole  based  on  mass  spectrometer  measurements 
given  by  Honig^  and  recommended  by  Stull  and  Sinke.  25  Searcy26  gives 
a  value  of  112  Kcalorie/mole  and  Huggins27  gives  100  Kcalorie/mole  for 
the  heat  of  sublimation.  The  results  of  a  very  recent  study  of  the  sub¬ 
limation  of  sic  by  Drowart,  deMaria,  and  Inghram2®  lead  to  a  value  of 
^  ^  Kcalorie  / mole  for  the  heat  of  sublimation,  in  good  agreement  with 
Searcy's  value  by  the  same  method. 

It  should  also  be  mentioned  that  a  very  recent  paper  by  Baird  and  Taylor2^ 
gives  a  heat  of  reaction  to  form  Sic  from  liquid  silicon  of  -19.  250  Kcalorie/ 
mole  as  compared  with  earlier  values  of  -23,  300  Kcalorie  /mole.  This 
leads  to  a  heat  of  sublimation  for  silicon  of  116  or  117  Kcalorie/mole 
from  the  results  of  Searcy  and  Drowart  et  al's  work. 

However,  these  high  values  (112  Kcalorie/mole  and  greater,  based  on 
Sic  dissociation  data)  should  not  necessarily  be  accepted  in  preference 
to  Honig's  value  which  is  based  on  a  direct  measurement.  A  further 
complication  in  the  measurements  based  on  SiC  dissociation  is  indicated 
Baird  and  Taylor's  statement  that  SiC  becomes  unstable  relative  to 
liquid  silicon  and  graphite  near  2,  050#  C. 

For  the  dissociation  energy  of  SiO(g)  ,  it  is  found  that  Barrow  and 
Rowlinson-50  recommend  185  ±  6  Kcalorie/mole.  This  agrees  with 
Gaydon's  value31  of  1  84  ±  23  Kcalorie/mole.  but  is  higher  than  the  value 
given  by  Herzberg  ^(171  Kcalorie/mole).  By  calculation  from  the  values 
of  105  Kcalorie/mole  for  the  sublimation  to  monatomic  silicon,  117.9/2 
Kcalorie  for  the  dissociation  of  one -half  mole  of  oxygen  from  Stull  and 
Sinke  3 and  a  heat  of  dissociation  of  185  Kcalorie/mole  for  SiO  ,  a 
value  of  -21,  100  calorie /mole  is  obtained  for  the  heat  of  formatfon  of 

^(j)  at  ab*olute  zero.  Other  possible  combinations  are  shown  in 
Table  II -A. 
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TABLE  H-A 


POSSIBLE  HEATS  OF  FORMATION  OF  SiO(g) 
FROM  SPECTROSCOPIC  DATA 
Heat  of  formation  of  Si°(f)  at  0*  K 

^n°f  ■  Lsi  +  «e>o2  “  DSiO 


Lo- 

Si 

Xl>02 

DSiO 

AH°r 

k&O) 

88  (ref 

19) 

58.  9  (ref  25) 

169  (ref  19) 

-22 

88  (ref 

19) 

58.  9  (ref  25) 

185  (ref  30) 

-38.  1 

100  (ref 

27) 

58.  9  (ref  25) 

169  (ref  19) 

-10.  1 

100  (ref 

27) 

58.  9  (ref  25) 

185  (ref  30) 

-26.  1 

105  (ref 

24) 

58.9  (ref  25) 

169  (ref  19) 

-  5.  1 

*105  (ref 

24) 

58.9  (ref  25) 

185  (ref  30) 

-21.  1 

112  (ref 

26) 

58.9  (ref  25) 

169  (ref  19) 

+  1.9 

112  (ref 

26) 

58.9  (ref  25) 

185  (ref  30) 

-14.  1 

*  Indicates  a  combination  of  data  which  agree  with  thermochemical  data. 

The  data  in  this  table  shows  that  a  wide  range  of  possible  values  for  the 
heat  of  formation  of  siO  could  be  estimated  depending  on  the  choice  of 
the  heat  of  sublimation  of  silicon  and  the  heat  of  dissociation  of  SiO  • 

In  view  of  these  uncertainties  it  would  appear  that  the  thermochemical 
data  reviewed  by  Brewer  and  Edwards®»^are  probably  more  reliable. 
Accordingly,  we  use  their  value  ofAH°0  -  -21,695  calorie  /mole  orAH^g  = 

—  21,411  ♦  574  calorie /mole. 

It  should  be  noted  that  the  data  of  Table  II -A  indicates  that  larger  heats 
of  sublimation  of  silicon  require  larger  heats  of  dissociation  of  SiO  in 
order  to  give  heat  of  formation  values  consistent  with  the  Brewer  value 
chosen  above.  Thus,  *'  a  value  of  117  Kcalorie  /mole  for  silicon  sublima¬ 
tion  were  to  become  accepted,  a  heat  of  dissociation  for  SiO  of  197 
Kcalorie  /mole  would  be  required  to  correspond  to  a  heat  of  formation  for 
SiO  of  -21.  5  Kcalorie /mole.  Such  a  value  for  the  heat  of  dissociation 


would  be  within  the  range  estimated  by  Gaydon  (184  ±  23  Kcalorie/mole.  ) 
At  this  time,  it  therefore  appears  that  thermochemical  data  are  to  be 
preferred  to  the  spectroscopic  data,  for  the  heat  of  formation  of  sio.  v 

(f)  * 

Having  now  selected  our  source  data  for  the  heat  of  formation  of  siOr  * 
at  lower  temperatures,  it  is  necessary  to  extrapolate  these  data  to  mgh 
temperatures  in  the  range  2,  000*  to  3,  000*  K.  To  make  such  an  extra¬ 
polation,  it  is  necessary  to  obtain  heat  capacity  data  for  SiO^  in  this 
temperature  range, 


To  obtain  thermodynamic  data  for  SiO(f)  ,  it  was  necessary  to  make 
calculations  with  spectroscopic  data.  The  methods  followed  are  similar 
to  those  described  by  Mayer  and  Mayer.  *3  The  molecular  weight  was 
taken  to  be  44.  09.  Other  spectroscopic  data  were  taken  from  Herzberg.  32 
No  electronic  levels  other  than  the  singlet  ground  state  were  of  import- 
***  these  calculations.  The  next  highest  electronic  state  has  a  wave 

n2Sjjrj?er  of  42»  835  «H-1,  corresponding  to  -^p-  =  20.6  at  3,  000#  K;  and 

e  makes  a  very  small  contribution  to  the  partition  function.  For 

the  one  vibrational  frequency  possible  in  this  diatomic  molecule,  the 
characteristic  frequency  (<u  -  <ue  -  2a»e*e  )  is  1,230  err.-1  .  For  the  rota¬ 
tional  contributions,  the  moment  of  inertia  was  taken  to  bel-38.6xl0~40 
I® cm2.  The  symmetry  number  is  one  since  the  atoms  of  the  molecule  are 
dissimilar.  Anharmonicity  corrections  to  the  entropy  were  made  at  the 
higher  temperatures  but  were  still  quite  small. 


The  results  of  these  calculations  are  tabulated  in  Table  II-B.  For  com¬ 
parison  purposes  the  results  of  Brewer  and  Edwards8*19at  the  lower 
temperatures  have  been  included  and  are  found  to  agree  with  our  data. 

A  further  approximate  check  on  our  data  is  indicated  by  the  fact  that  it 
was  found  that  the  heat  capacity  term  due  to  vibration  at  2,  000*  K  was 
almost  fully  excited.  Moreover,  if  for  the  diatomic  gas  we  assume  that 
transl&tioni  rotation,  and  vibration  are  contributing  fully  to  the  heat 
capacity  at  constant  pressure,  we  get  C.-9/2R  .  Therefore 

x  om  r 


h'OOO  -  ^.ooo 


/ooo 

2,000 


Cp<fr  ^ 
T 


3*62  c.a 


From  our  results  in  Table  II- A  the  entropy  increment  is  3.63  e.u.  Calcu¬ 
lations  with  the  enthalpy  data  also  lead  to  _  1*3000^.876  calorie/mole, 

giving  an  average  0  over  this  temperature  range  of  8, 876  calorie  /mole 
degree.  Thus  the  data  appear  to  be  self-consistent. 


j 


With  the  data  from  Table  II- B  and  additional  data  for  elemental  silicon 
and  oxygen  from  Stull  and  Sinke2?  it  is  now  possible  to  calculate  the  heat 
of  formation  of  gaseous  siO  in  the  temperature  range  of  2,  000*  - 
3,  000»  K.  The  results  of  such  a  calculation  are  summarized  in  Table  m. 

The  calculation  was  based  on  the  relation 


( 

♦ 

l 

I 

I 

I 

I 


or 


o 


-  ah. 


♦  <Ht  -  Holsio^  -  <Ht  -  H^si  -  K(Ht  -  1^)0 


2(*> 


In  order  to  make  an  estimate  of  the  uncertainty  in  this  calculation,  it 
was  necessary  to  take  into  account  the  accuracy  of  all  the  data  used  in 
making  it.  Since  Stull  and  Sinke25had  not  stated  the  accuracy  of  their 
ta,  we  made  our  own  estimate  of  their  accuracy  as  indicated  in  Table 
III.  In  order  to  make  a  partial  check  on  our  computation  of  the  heat  of 
formation  of  SiO(f)  .  we  have  tabulated  the  data  of  Coughlin9  and 
Humphrey  et  al  in  the  second  row  of  the  table.  Their  data  for  298*  K 
and  2000*  K  are  in  good  agreement  with  ours. 


3*  Heat  of  the  Overall  Reaction 

In  the  previous  sections,  the  heats  of  formation  of  siO,,  .  and  siO 
in  the  range  of  2.  000*  K  -  3.  000*  K  were  evaluated.  T Zoie  reeult.  c£i 
be  readily  used  to  calculate  the  heat  of  the  overall  reaction  for  the 
vaporization  of  liquid  siC>2  by  mean*  of  the  relation  AH,..  AH,  [siO^J-  AH,  [s»2(1)] 

^5T00-1K.IVandrVooT1K.riZed  ^  r"Ult'  °£  *UCh  **  2.  °°°* 


J 
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TABLE  II-B 


THERMODYNAMIC  FUNCTIONS  FOR  Si0(^  FROM  SPECTROSCOPIC  DATA 


ALL  UNITS  ARE 

ft' T  “  «•  'N 

298*  K 

\  T  S 

Translation 

-32.  30 

Electronic 

0.  00 

Vibration 

0.  00 

Rotation 

-11. 27 

ft°  -  H*\ 

l  T  °) 

\  T  /“ 

-43. 57 

("-*) 

\  T  / 

T  ranslation 

4.  965 

Electronic 

0.  000 

Vibration 

.  00 

Rotation 

1.  987 

/H  - 

6.  955 

50.  52 

net  “« 

Anharmonicity 

Corrections 

0.  00 

ST  (our  calculation) 

50.  52 

S  (from  Brewer  and 

T  Edwards®*  19) 

50.  53 

CALORIE/ MOLE -DEGREE  K 


2, 000#  K 

2,  500* K 

3,  000*  K 

-41.76 

-42. 86 

-43.76 

0.00 

0.  00 

0.00 

-  1.05 

-  1.34 

-  1.61 

-15.05 

-15.49 

-15.85 

-57.86 

-59.69 

-61.22 

4.  965 
0.000 

1.23 

1.987 

4.965 

0.  00 

1.36 

1.987 

4.  965 
0.00 
1.46 
1.987 

8.  18 

8.  317 

8.412 

66.  04 

68.  01 

69.  63 

0.  00 

0.  03 

0.  04 

66.  04 

68.04 

69.  67 

66.  02 
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TABLE  IV 


HEAT  OF  THE  REACTION  Si02(i)— ►  SiO(g)  +1/2°  2(g) 


2, 000*  K 

2, 500*  K 

3,  000*  K 

AH* 

-36,  135  ±  1, 264 
or  (±  2,  764) 

-37, 527  ±  1, 394 
or  (±  2,894) 

-38,  954  db  1,  543 
or  (±  3,  043 

AH# 

'<»°2(1)> 

AH 

reaction 

-215,  200  ±  1.  000 

-212,  453  ±  1, 615 

-210, 467  ±  2,  230 

179.  065  ±  2,  264 
or  (±  3,690) 

174, 926  ±  2,  834 
or  (±  4,  260) 

171,  513  ±  3,  409 
or  (±  4,  835) 

Note  that  uncertainties  at  2,  500*  K  and  3,  000*  K  are  less  than  would  be 
calculated  by  simple  addition  of  the  component  uncertainties.  (See 
Table  V.  )  A  second  set  of  estimated  uncertainties  (indicated  in  paren¬ 
thesis)  are  1, 500  calorie  in  excess  of  those  not  shown  in  parenthesis. 
Thus.  Coughlin  ^  gave  an  uncertainty  of  2.000  calorie  for  AHf(Si0  )  at 

298*  K,  whereas  Brewer  and  Edwards  gave  574  calorie  even  though  they 
had  both  evaluated  the  same  original  data.  To  make  the  best  estimate 
of  the  uncertainties  in  the  heats  of  reaction  obtained  above,  the  sums 
involved  in  deriving  the  values  must  be  considered  in  detail.  It  is  shown 
by  the  following  algebraic  manipulations  that  some  of  the  uncertainties 
cancel  out.  The  remaining  terms  actually  used  in  calc  ulating  the  heat  of 
reaction,  can  be  added  without  regard  to  algebraic  sign  to  get  the  true 
uncertainty. 

For  the  reaction  SiOjQj— ►  SiO^  ♦  it  follows  that 


AHMrtion  -  AH,  +  AH 
reaction  q  j 


f(siCY 


-  O  ♦  AH' 


'<SiO(t)> 


♦  (HT  -  H^sio^  -  (Ht 


-  Ho>Si(1) 


-  ^Hr-H^o 


-  (HT  -  h2ooo>  ♦  (HT  -  H2000)  ♦  (ht  -  h2000>, 


SiO- 


Si 


(1) 


*<•> 
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or 


^reaction  “  ^pf  “  A*1 2 000 f  *  WT  -  Hj,  )  -  (HT  -  H^) 

'Si°(«)>  *Si°2,  J  (Si°(*)>  W(SK>,_.) 


(  1) 


2(1) 


“  (H2000  ”  ^g.  ♦  l^(HT  -  H2000>0  -  ,/4(H2000  -  H,,) 

(1)  2(«)  2(*>  . 


Table  V  shows  the  magnitudes  of  the  uncertainties  of  the  quantities  in  the 
latter  equation.  Their  stuns  were  taken  as  the  true  uncertainty.  At 
3,  000*  K  it  is  seen  that  about  half  of  the  error  arises  from  the  thermo¬ 
chemical  data  for  the  heats  of  formation  of  the  reference  materials.  This 
may  be  somewhat  optimistic  inasmuch  as  the  other  determinations  of  the 
heat  of  formation  of  siO^  j  quoted  by  Brewer  and  Edwards^* ^indicate 
a  possible  spread  of  a  feVr  thousand  calories  per  mole.  Hence*  it  is 
possible  that  the  uncertainty  in  AH®  may  be  larger  than  the  574 

°Si0<«> 

calorie /mole  given.  If  this  uncertainty  is  2,000  calories  /male,  as 
chosen  by  Coughlin,  then  the  resultant  uncertainty  in  the  overall  heat  of 
reaction  at  3,  000*  K  is  about  5,  000  calories /mole.  The  2*000  value  for 
the  uncertainty  of  the  heat  of  formation  of  siO(  >  may  be  more  reason¬ 
able  than  the  574  value  since  other  experimenters'  values  have  had  a 
wider  spread  than  574  calorics.  The  other  large  uncertainty  contribu¬ 
tions  arise  from  the  estimation  of  heat  capacity  data  for  liquid  Si02 
and  liquid  silicon. 

For  the  previous  calculation  of  the  heat  of  the  reaction  siO  — ►SiO  +$40 

it  might  appear  preferable  to  use  heat  of  reaction  data  base£(£>n  this(f> 
reaction  directly.  However,  directly  determined  data  such  as  those  of 
Brewer  and  Mastick2  and  recalculated  by  Brewer  and  Edwards8,  *9  giveB  a 
greater  uncertainty  in  the  final  heat  of  reaction.  For  example,  these 
authors  give  AH°.  186,578  +  2,700  calorie/mole  for  the  reaction  SiO-,  . 

— ^  +  ^2(a)  •  a  check  on  the  earlier  calculations  and  in  orae?,#t 

to  show  that  the  Brewer  Mastick  work  is  more  uncertain,  it  is  possible 
to  calculate  A  H-j-  for  the  first  reaction  above  using  the  direct  data  of 
Brewer  and  Mastick  in  the  following  manner: 


AH. 


-  ah°c 


♦  (H, 


“  Ho^SiO 


(*> 


♦  *<HT  -  “obxj  -  <H 


“  Ho^SiO 


2(1)  referred  to 

cnetob&lite  at  0°K 
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TABLE  V 


SUMMARY  OF  COMPONENT  UNCERTAINTIES  OF  THE  HEAT  OF  REACTION 

All  units  are  calorie /mole. 


Uncertainties 

2, 000  K 

sH,) 

±  574  (2.  000)^ 

^iooot 

*°2U) 

±  1.000 

<Ht  -  *  w> 

±  150 

(HT  -  H2000>  iQ 

Sl°2(l) 

0 

(H2000  -  . 

*(1) 

±  460 

*<«T  ~  Ho) 

°2(«) 

H- 

00 

o 

±  2.,  264(3,  690)++ 

2, 500  K 

3, 000  K 

±  574  (2,  000)^ 

±  574  (2,  000)^ 

±  1, 000 

±  1,000 

±  200 

±  250 

±  500 

±  1,000 

±  460 

±  460 

±  100 

±  125 

±2,  834(4,  260)** 

±  3,  409(4,  835)** 

♦Combined  uncertainty  of  last  two  terms  of  equation  used. 


♦♦Coughlin's  estimate  of  uncertainty  in  AHf/cjoi 
is  greater  than  that  of  Brewer  and  Edwards. 


shown  in  parenthesis. 
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Value*  of  (Ht  _  k^o 

obtained  and  listed  in  t 
to  cristobalite  at  0*  K, 


and  <ht  “  h°)°2 (  »  have  already  been 

able  III*  The  heat  content  of  liquid  Si 02  relative 
can  be-obtained  from  the  relation 


(HT  “  ^SiO^j)  "  <«T  -  H2ooo)si02(1)  ♦  AH 


‘fusion  of  * 
cristobalite 


(H 


2000  “  H»8)criJH»8  ~  Ho>  . 

crtmt.  cnst. 


which  can  be  evaluated  as  follows: 


(HT  “  H2000>  . 


SiO 


2(1) 


^fusion 


(H2000  “  H298) 

cnst. 


-  21,660  ±  1,000  calorie /mole  (from  Table  I.) 

=  2,500  ±  500  calorie/mole  (see  earlier  dis¬ 
cussion) 

=  28,  120  ±  600  calorie /mole  (K.  K.  Kelley13) 


(Hjoo  —  H  ) 

298  °  crist. 


=  1 , 667  ± 


(ht  -  «.)  . 


SiOjf.v  =  53>  947  ± 

‘'referred  to 
cnstobslite  at  0°K 


50  calorie /mole  (Edwards19) 
2,  150  calorie /mole 


Substituting  the  latter  into  the  expre.aion  for  AH  .  we  find  that  AH 
*  170.621  ±  5,225.  T 


5000 


!  **  .ft  J ’ Vlooti0n  agree*  Wel1  wlth  our  previous  value  of 
of  tI'mb  IV  h.7  ’  513  ,,3’  400  Calorle/mole-  However,  the  earlier  data 

cullttons.  <,maUer  “"“rtainty  and  are  used  for  subsequent  cal- 

Total  Heat  of  Decomposition 

The  total  heat,  which  can  be  absorbed  bv  fused  . 1 1 , 

at  3,000*  K  bv  the  reaction  _  V  fU“*d  *lllCa  when  “  decompose. 


at  3,  000*  K  by  the  reaction  s^q 
as  follows : 


(g)  +  ^°2(g)  *  ia  now  computed 


(”2000  ~  ”298><yo 

(glsss) 


=  29,010  calorie/mole  (Kelley13) 
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*  H3000  “  H2000>SiO 


2  (liq  or  flats) 


=  21,660  calorie/mole  (Table  I) 


AlljOOO  *»»  sio2 — ►  SiO(f)  ♦  xo2 

Total  heat  of  decomposing  SiO, 
at  3,  000*  K  starting  from  ***•■•> 
298*  K. 

The  equivalent  value  in  other 
units  at  4.  940*  F  starting  from 
77*  F  is 


B.  ENTROPY  CHANGE 


=  171,513  calorie/mole  (Table  IV) 


=  222,  183  calorie /mole 


=  6,  650  Btu/lb 


Stull  and  Sinke25give  the  absolute  entropy  of  02(  *  .  in  Table  H  are  given 

the  entropy  data  for  SiO(f)  .  Kelley*  gives  We  entropy  of  at 

298*  K  as  s^g  =  H.2  ±  0.4  eu.  For  higher  temperatures  Kelley13gives 


S2000  “  s298 


=  29.  87  eu. 


entropy  data  for  sio 


ThU*  S2000 


=  41.  07  eu.  In  order  to  obtain 


.  .  J  *2(1)  above  2, 000*  K,  we  have  had  to  use  our  estimated 

heat  capacity  data  (c  =21.66  calorie/mole  degree).  The  results  of  these 
calculations  are  sho*n  in  Table  VI. 

C.  FREE  ENERGY  CHANGE 

Having  obtained  the  heat  of  reaction  and  the  entropy  change  in  the  previous 

sections,  it  is  now  possible  to  calculate  the  free  energy  change  from  the  re¬ 
lation 


AF  -  AH  -  T  A  S 


The  results  of  this  calculation  are  shown  in  Table  VII. 

For  AH  »  the  uncertainties  represent  maximum  estimates.  In  estimating  the 

uncertainties  in  AF  from  the  data  in  Table  IV,  the  uncertainties  of  AH  and 
TAs  do  not  add  algebraically.  Actually.  *’  ->se  uncertainties  cancel  in  part 
by  the  use  of  consistent  heat  capacity  data  ....  discussed  in  the  following  para- 
graphs  so  that  the  value  given  at  3,  000*  K  is  probably  good  to  about  *  4.  000 

calories /mole,  depending  again  on  the  uncertainty  of  the  heat  of  formation 
of  SiO(([)  . 
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TABLE  VI 


ENTROPY  CHANGE  FOR  THE  REACTION  Si02(1)— siO(g)  ♦  X02(g) 


Temperature  in  Degrees  Kelvin 


2,  000*  K 

2, 500*  K 

3,  000*  K 

Ref. 

'AS°2(t> 

(e.  u.  ) 

32.  1 1  db  .02 

33.  13  ±  .02 

33.  99  ±  .03 

2 

Ssi°(,) 

(e.  u.  ) 

66.  04  ±  0.  1 

68.  04  ±  .1 

69.  67  ±  .1 

Table  II 

Ssi°2(l) 

(e.  u.  ) 

41. 07  ±  0.  6 

46.  13  ±  0.  8 

50.  03  ±  1.0 

1.  9  at  2,  000*  K 
calculated 
using  Cp  =  21. 66 
calorie/mole  #  K 
for  T  > 

2, 000*  K 

ASj* 

(e.  u.  ) 

57.  08  ±  0.  7 

55.  04  ±  0.  9 

53.  63  ±  1.1 

TAS  calorie  /mole 

114,  160  ±  1, 400 

137,  600  ±2,  250 

160,  890  ±  3,  300 

TABLE  VII 

FREE  ENERGY  CHANGE  FOR  THE  REACTION  Si°2(1)  — ♦  X°2(g) 

Z,  000*  K  2, 500*  K  3, 000* K 


AH 

(calorie  /mole) 

179.  065  ± 

3,  700 

174, 926  ± 

4,  300 

171, 513  ± 

4,  800 

-  TAS 

(calorie  /mole) 

-114,  160  ± 

1, 400 

-137,  600  dt 

2,  250 

-160, 890  ± 

3,  300 

af 

(calorie  /mole) 

64, 905  ± 

3,  000 

37,  326  ± 

3,  500 

10,  623  ± 

4,  000 
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Th®  variation  of  the  above  free  energy  data  with  temperature  can  be  repre  - 
sented  with  sufficient  accuracy  by  a  straight  line  such  as 

AFT  -  173*000  —  34*2  T  (caloric/ mole) 


which  is  valid  for  the  temperature  range  of  2.  000*  K  to  3,  000*  K  and  accurate 
to  ±  4,000  calorie/mole  at  3,  000*  K,  The  values  from  this  equation  fit  the 
calculated  points  of  Table  VII  sufficiently  well  that  it  can  be  used  for  the  sake 
of  simplicity  in  our  applications.  This  linear  equation  can  also  be  used  to 
extrapolate  slightly  beyond  the  temperature  range  of  the  data.  (2.000  - 
3,  000*  K) 

In  assessing  the  accuracy  of  the  free  energy  expression  which  we  have  ob¬ 
tained,  we  may  take  as  a  second  trial  a  heat  capacity  for  liquid  siO  of 

20.  66  calorie /mole  •  K  instead  of  the  21. 66  calorie /mole  •  K  as  in  34V  earlier 
calculations.  With  such  a  value  of  heat  capacity,  the  heat  of  formation  of 
Si°2(l)  at  3*  00°*  K  would  be  AHjoqo*  =  -211, 467  (instead  of  -210,  467) 
calorie  /mole.  In  Table  IV  use  of  this  value  of  AH.  will  raise  the 

heat  of  the  reaction  to  AH  [si02(1)—  SiO(g)  ♦  «  °2(, J’  2<I>  =  172.  513 

calorie/mole  (instead  of  171. 513).  Next,  from  Able  VI  it  can  be  seen  that 
the  entropy  of  would  be  lowered  by  the  use  of  the  smaller  Cp  value 

by  an  amount  equal  to  AS^O  n  =  0.  41  e.  u.  This  would  make  ^ 

SiOj-  ,  000° If 

=  49.62  e.  u.  (instead  of  50.03  e.  u.  ).  Consequently,  the  As  of  reaction 'at 
3.  000*  K  would  be  increased  to  54.  04  e.  u.  (instead  of  53.  63  e.  u.  ),  and 
(T As) 3000° k  =  162,  120  calorie  (instead  of  160,890  calorie).  The 

corresponding  free  energy  value  is  AF3f0()0oK  =  10,  393  calorie /mole 

instead  of  the  10,623  calorie  /mole  obtained  with  q,  =  21.66  calorie /mole*  K 
or  a  decrease  below  the  previous  value  of  2  30  calorie /mole.  The  lower  c„ 
f°r  SiOj  indicates  a  slightly  increased  vapor  pressure. 

It  is  seen  from  the  above  discussion  that  the  free  energy  does  not  change  too 
rapidly  with  changes  in  the  value  of  the  heat  capacity  of  liquid  silica  although 
its  components,  the  heat  of  reaction,  and  the  entropy  do.  This  can  be  shown 
by  the  following  more  general  treatment. 

To  estimate  the  accuracy  of  our  free  energy  functions,  we  need  to  reconsider 

in  detail  how  they  were  calculated.  The  discussion  is  based  on  the  fundamental 
definition, 

AF  -  AH  —  TAS  . 
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In  considering  the  reaction 


si°2(l>  — ►si°(g)  ♦  ^°2(*) 


(7) 


we  have  noted  that  at  2,000*  K  the  free  energy  change  is  accurate  to  about 
±  3,  000  calories  /mole.  This  is  the  uncertainty  estimated  from  Coughlin's 
table.  Of  this  uncertainty  about  2.  000  calorie/mole  is  associated  with  con- 
tnbutions  from  the  formation  of  SiO(ff)  and  1.000  calorie/mole  is  associated 
With  the  formation  of  liquid  Si02  .  To  make  a  rough  estimate  of  the  errors 

to  be  expected  at  higher  temperatures,  we  can  modify  the  fundamental  equa¬ 
tion  to  obtain 


AF  -  AH 


2000 


2000 


i 

J  AC,,  dT  - 
2000 


,<1T  -  T  ^ 

T 

As2000  ♦  /  ACp  dT  | 
2000  T  / 

2000  ~  (T 

—  2000)  AS 2000 

T 

/ 

2000 

ACpdT 

T 

If  we  now  designate  the  first  two  terms  of  this  last  equation  as  Af  an 

replace  our  variable  of  integration  T  by  r  the  dummy  variable  (to^void 
confusion  with  the  T  used  as  the  limit  of  our  integration)  we  obtain 


AF  “  A  F 2000  ~  (T— 2000)  ASpOOO 


r 


dr 
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Af2000  “  (T  —  2000)  ASjoqo  ♦  T 


( AQ»  _  AQ>^  ^ 
T  r 


The  uncertainty  in  Af2#00 0  is  ±  3,  000  calorie /mole  according  to  our  earlier 
estimate.  Since  the  value  of  AS2,000  =  57.  08  ±  0.  7  e.  u.  .  the  uncertainty 

in  the  second  term  at  3,  000*  K  is  ±  700  calorie /mole.  The  last  integral 
would  also  contribute  to  the  overall  uncertainty  if  ACp  for  the  reaction  were 
in  error.  For  the  present  case,  we  may  get  a  rough  order  of  magnitude  of 
this  uncertainty  if  we  note  that 


<=P  for  %02  ^  4.  5  calorie/mole  degree,  =  9.  0  calorie /mole  degree. 

and  Cp^  =  21.  66  calorie/mole  degree.  Hence.  Aq>  'T  -8.0  calorie/ 

mole  degree.  The  most  likely  source  of  error  in  Ar^  is  .If 

error  were  1  calorie/mole  degree  so  that  actually  Cp«yo  ^  =  20.66 

and  ACp  =  -7.  q,  then  the  total  error  would  amount  to  an  additional  220 

calorie  /mole.  Thus  at  3,  000*  K  the  uncertainties  in  AF  add  up  to  3,  920  or 
4,  000  calorie  /mole. 


As  a  further  check  on  the  free  energy  equation  which  we  have  obtained,  we 
can  compye  values  of  the  free  energy  at  2,  000*  K  with  those  tabulated  by 
Coughlin.  He  gives  AF®^  =  -59.  200  calorie/mole  and  Af^qoo 

Si°(g)  *  s‘°2(l) 

”  -123,900  calorie/mole  from  which  we  get  for  the  reaction  siO  _ ►  siO,  x 

2(d  Cg) 

♦  l/t°2(g)  t*ie  *ree  energy  change 

AF2tooO  •  64,700  £  3*000  calorie/ mole 


This  compares  favorably  with  our  value  of  64,600  calorie/mole  based  on  the 
empirical  formula. 

In  figure  1  are  plotted  the  free  energy  data  obtained  by  the  present  calcula¬ 
tions.  The  empirical  equation  has  been  extrapolated  above  3,000*  K.  This 
can  be  safely  done  over  a  narrow  temperature  range  and  will  be  a  fair  approxi 
mati on  at  higher  temperatures.  From  the  expression 
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AF  -  AH  -  TAS  , 


with  a  reference  temperature  of  3,  000*  K,  we  find 


T  T 
AFr  "  AH3.000  +  J  AC^dT  -  TAS3t000  -  T  f 

3.000  3,000 


ACpdT 

T 


where 


i.  e.  , 


ACp  is  equal  to  the  change  in  heat 


capacity  of  the  reaction  considered. 


Acp  ■  ‘AcPo2  *  CPSiO  -  SsiO 


2(1) 


For  .mall  temperature  intervals,  the  second  and  fourth  term,  in  the  free 
energy  equation  will  cancel,  thus  justifying  the  extrapolation. 


AFT  -  AH^ooo  -  TAS3  00O 


For  temperature,  below  2,  000*  K.  where  solid  SiO,  is  tile  stable  phase 
the  data  of  Coughlin  for  cristobalite  have  also  been  plotted  in  figure  1.  'in 
this  case  hi.  tabulated  free  energies  of  formation  for'  SiO,  .andV  can 
be  subtracted  and  used  directly. 


It  is  also  possible  to  develop 
diate  range  below  2,  000*  K, 
the  reaction 


an  empirical  free  energy  equation  for  the 
by  the  use  of  entropy  data.  For  example. 


imme  - 
for 


^*^2 (cristobalite)  +  ^®2(g) 


we  get 


^cristobalite;  2,000o  k  "  S298°K  +  (s  2,000°K  ~  S 


>298°K 


:j  - 


•04  e.u.  13,14 


H  Sq 

2;2,000©k 

SSiO  2,000 °K 


32.11  e.u. 
66.O4  e.u. 
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and 


As2,000°K  "  SsiO  ♦  ^So2  “  %i02  "  58.Jle.u. 


Using  this  value  of  the  entropy  change  at  2,  000*  K,  we  can  now  fit  a  curve 
through  the  2,  000*  K  point,  since  AF2f000  =  64,600  calorie /mole  from  our 

previous  formula.  Hence, 

AF°  «  180,800  -  58.11  T(c*lorie/mole  for  1,800  <  T  <  2.000°K) 


The  above  relation  gives  A F°i  a00  =  76,200  calorie /mole  and  tot°1900 

=  70,  390  calorie /mole.  These  values  can  be  compared  with  values  calculated 
from  Coughlin's^  data 

A F°h  000  “  76,450  calorie/mole 
^  1,900  —  70,650  calorie/mole. 

This  good  agreement  increases  our  confidence  in  the  data  and  calculations. 

D.  EQUILIBRIUM  CONSTANT 

For  the  evaluation  of  the  equilibrium  constant  corresponding  to  the  reaction 
that  is  the  subject  of  Section  II,  the  usual  relation  AF°  =  -RT  in  K  can  be 
modified  to  give 

-  AF  AH  AS 

lotinK  m  -  ■  ■  *  —■  • " 

1  2.303RT  4-57T  4.57 


For  the  present  case,  this  becomes 

,  „  -37,900 

Io«l0K  "  — - -  +  11.85 


This  equilibrium  constant  has  been  plotted  in  figure  2  as  a  function  of  1/T.  It 
can  be  seen  that  at  higher  temperatures,  the  equilibrium  constant  becomes 
larger. 
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E.  EQUILIBRIUM  VAPOR  PRESSURE 

In  the  range  2,  000  -  3,  000*  K  we  have  ihown  that  the  standard  free  energy 
o  e  reaction  SiOjQ)  can  be  represented  as  a  function 

of  temperature  by 

AF|  -  173,000  -  54.2T. 

This  is  related  to  the  equilibrium  constant,  K  by  the  conventional  definition 

AF  ^  --RTlnK  -.RTln  -l0<S>  *°2(g) 

*Sl°2(  1) 


where  =  activity  of  component  (or  the  fugacity  in  case  of  gaseous  compo 
nents).  For  liquid  SiC>2(1)  ,  the  usual  assumption  is  made  that  the  activity 

is  unity,  i.e.  ,  any  solution  effects  of  SiO(g)  #  or  other  equilibrium 

species  in  S>02  are  considered  negligible. 


For  the  gaseous  components,  it  is  reasonable  to  replace  the  activities  by  partial 
pressures  of  the  species  being  considered  as  long  as  the  gases  do  not  deviate 
markedly  from  the  perfect  gas  laws.  This  assumption  is  valid  at  low -pressures 

n"r-hjfth’t*niperatUre8*  F°r  OXy8en  (fugacity  =  48  atm  at  pressure  of  50  atm  and 
0  C.  critical  temperature.  Tc  =  154.  3  •  K,  critical  pres  sure.  Pc  =49.  7atm)at 


50  atm  and  2.  000*  K  the  reduced  temperature  is  ^  Z 

R  *r 
Ac 

the  reduced  pressure  Pf  -  p/pc  =  0.  995.  For  this  case34 

and  the  fugacity  is  essentially  equal  to  the  pressure. 


=  12.  95,  and 
=  1.0  + 


For  SiO(^)  ,  whose  condensed  properties  are  practically  unknown,  we 
must  use  some  very  rough  approximations  to  estimate  pc  and  Tc 
Hirschfelder  £t  al  list  some  rough  empirical  rules  which  can  be  used  in 
this  case.  Using  the  condition  that  the  melting  temperature  is  2/5 
and  the  fact  that  EdwardsiV  states  that  the  melting  point  of  sio  must  be 
greater  than  1,975*  K,  we  calculate  Tc  ^5,000*  K.  To  estimate  the  molar 
volume  at  the  crxtical  point,  we  note  that  V*®1  =  0.321  .  According 

t0  ^dwards  the  density  of  solid  SiO  i8  about  2.  15  gm/cm?  (20.  5 

cm  /mole).  Making  the  assumption  that  the  molar  volume  of  siO  at  its 
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melting  point  is  equal  to  the  molar  volume  (20.  5  cm3/mole)  given  by  Edwards, 
we  find  a  molar  volume  at  the  critical  point  of  ~  =  64.  0  cm3/mole. 

Finally,  from  the  fact  that  many  mplecules  have*a  compressibility  factor  of 
0.  292  ( Hirschfelder  et  alf2  ,  we  find  that  pc  =  1,870  atm.  At  50 

atm.  and  2,000*  K,  p^  =  0.0267,  and  T^  =  0.4.  Again,  estimating  from 
the  Hougen  and  Watson  chart,  it  is  found  that  f/p  »  0.  95.  In  view  of  the 

uncertainties  involved  in  the  above  calculations,  the  results  must  be  con¬ 
sidered  very  approximate.  However,  they  do  indicate  that  in  the  pressure 
range  of  interest  (1-50  atm)  and  from  2,  000  -  3,  000*  K,  the  gas  will  behave 
ideally.  For  further  considerations,  pressures  will  be  used  instead  of  fugaci- 
ties. 

Using  the  simplifying  assumptions  discussed  above,  we  get  for  the  final 
formulas 


RT  In 


A 

PSIO  P02 


RTlnKeq  , 


or 


PsiO  POz 


-AF°  ->173,000  Mil 
RT  „  e  RT  e  R 


Under  neutral  conditions, 
and 


PsiO  ■  2Pq2 


PSiO 


/  3160* 
18. <1  (l  -  — — ) 

-  «  x  T  *  . 


The  vapor  pressure  data  for  this  equation  are  presented  as  the  upper  curve 
in  figure  3.  It  gives  a  partial  pressure  of  0.667  atm  for  siO(  at  3,  085*  K 

(Boiling  point  under  neutral  conditions,  neglecting  other  reactions).  It  should 
also  be  noted  that 


P*ot*l  *  1,3  pSiO  "  13 


18.41  (l- 
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The  vapor  pressure  of  SiO^j  in  equilibrium  with  condensed  Si02  is  plotted 
for  a  large  temperature  range  in  figure  4.  The  principal  curve  in  this  figure 
is  based  on  the  single  reaction  given  and  neglects  all  other  equilibria. 

For  the  case  of  Si02  decomposing  in  the  presence  of  air, 

PN2  “  4(p02  -  PsiO) 


and 


Ptotal  “  PN2 
From  the  above 

we  now  find 


♦  PsiO  “  constant  -  5po2  -  PSjO  * 
relation  for  psiQ 


Psi° 


'A 


e28.0  (1  -  3110) 
T 


These  data  have  been  plotted  in  figures  3  and  4  for  the  two  cases,  p  =  1  atm, 
and  p  =  1 0  atm. 

It  can  be  seen  from  the  above  equations  that  the  presence  of  the  oxygen  sup¬ 
presses  the  SiO  gas  formation.  It  should  also  be  emphasized  that  the  data 
here  presented  assumes  only  the  reaction  SiO(g)  +  l^°2(g)  to  be  important. 

Dissociation  of  oxygen  has  been  neglected  here  but  is  included  in  later  dis¬ 
cussions.  Nitrogen  is  considered  to  be  an  inert  gas. 

For  temperatures  below  2,  000*  K  the  relation  for  the  vapor  pressure  p  .n 
must  be  modified.  Using  the  relation  AF°  =  180,800  -  58.  11  T  for  * 

(1 , 800°  K<T <2,  000c  K) 

we  get 


19.72  (1  —  3  973) 
PsiO  *  *  T 
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A  curve  for  this  function  has  been  plotted  in  figure  5  over  the  temperature 
range  (1,800  -  2,  000*  K).  For  comparison  purposes  the  available  experi¬ 
mental  data  have  also  been  plotted*.  The  effusion  cell  measurements  of 
Brewer  and  Mastick  are  seen  to  agree  well  with  the  c  lculated  data.  The 
mass  spectrometric  measurements  of  Porter  et  al  3  are  h*gh  by  a  factor 
of  about  6.  Porter  et  al  indicated  in  their  paper  and  ,.n  j  ri,ate  communi¬ 
cations  that  their  absolute  values  of  pressures  could  be  in  error  by  a  factor 
of  2  or  3.  A  point  of  possible  uncertainty  in  the  work  of  Porter  jet  al  is  the 
fact  that  they  used  some  older  vapor  pressure  data  from  Brewer.  They 
had  calibrated  their  equipment  by  assuming  that  pA107  =  1  x  10~6  atm  at 

1.24S*  K#  and  that  pAu  =  2  x  10"5  atm  at  1. 750*  K.  however,  more  recent 
work  by  Honig 1  gives  values  of  1.  85  x  10“6  atm  for  a*107,  and  7.  9  x  10"5 
atm  for  Au.  All  observed  pressures  of  the  various  species  are  then  increased 
for  the  appropriate  temperature  range  by  the  ratios  1.  85/1  and  7.  9/2.  0.  It 
is  not  completely  clear  that  the  pressures  should  be  increased  by  this  much 

since  Porter  et  al 3  state  that  they  checked  their  calibration  by  the  vaporiza¬ 
tion  of  a  weighed  amount  of  Ag  or  Au. 


1' m TS?”  ,  data  fur,her  with  present  result.,  we  have  shown  in 
Table  VII-A  free  energy  data  calculated  from  their  observed  vapor  pressures. 
It  is  seen  that  the  result,  based  on  their  tabulated  pressures  are  in  good 
agreement  with  our  free  energy  data.  It  should  also  be  noted  that  their  reac¬ 
tion  does  not  correspond  to  stoichiometric  conditions.  They  apparently  did 
not  have  neutral  conditions.  The  free  energy  data,  ba.ed  on  recalculated 
vapor  pressures,  are  in  poorer  agreement  with  our  calculated  free  energy 
curve  given  in  figure  1.  Further  evidence  that  Porter's  tabulated  data  are 
probably  preferable  will  be  produced  when  the  dissociation  of  oxygen  is  con- 
tidered  in  a  following  section. 

In  the  literature  given  by  Inuzuka  and  Ageha38  are  some  additional  vapor 
pressure  data  for  Si02  .  They  report  a  vapor  pressure  ofl.lx  lO"7  atm 

°r  ^  k  ^  20°*  °  0  *  4?3*  They  used  a  "Mo  and  W  wire”  to  heat  a  fused 

Si02  tube  and  measured  the  vapor  pressure  by  a  Langmiur  method.  At  this 

1T?0-10 th*rmo^hemical  calculation,  give  a  total  vapor  pressure  of 
V*  10  .  fig.  4),  a  discrepancy  by  a  factor  of  180.  However,  it  is 

***•  f r«ument.  that  ^ey  were  actually  studying  some  other 

reaction,  possibly  the  following: 


1/3  M»(e)  *  SiOJ(c)  — ►SiO(j)  ♦  1/3  MoO,(j)  . 
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TABLE  VII-A 


COMPARISON  OF  FREE  ENERGIES  BASED  ON  DATA  OF  PORTER  et  al  3 


1, 800*  K 

^SiO(ntm) 

Po  „(•«») 

2 

^^(cslcd) 

(calorie  /mole) 

^^(our  curve) 
(calorie  /mole) 

76, 000 

Porter  data, 

as  is 

6.  8  x  10"6 

5.2  x  10“8 

72, 400 

Porter  data, 
corrected 

2.68  x  10-5 

2.  05  x  10"7 

65,  100 

1.  900*  K 

Porter  data, 
as  is. 

3.  4  x  10“5 

8.  2  x  10’8 

69, 595 

71,  000 

Porter  data, 
corrected 

1.  34  x  10“4 

32.  4  x  10'8 

61, 800 

than 


This  reaction  would  tend  to  give  a  higher  vapor  pressure  of  siO( 

expected  under  neutral  conditions.  The  presence  of  W  might  alto  reduce  siO  , 

but  it  would  not  be  expected  to  change  the  calculated  pressure  due  to  Mo  by  2 
more  than  a  small  factor.  y 

first  utin,thefrm°ChemiCal  da‘a  f°r  th*  reaction  with  Mo.  we  can 

first  us e  the  free  energy  data  of  Coughlin  for  the  reaction  to  give  liquid  MoO* 
instead  of  the  gaseous  product  q  Mo°  3’ 


l/3Mo(c)  +  Si02(c)  — +SiO(g)  +  1/3  Mo03(1) 
The  free  energy  of  this  reaction  is  then 


AF 


1473  -  A^siO  +  J/3  AF 


MoO, 


-  AF 


SiO- 


"  63.000  cal/ mole 


•Coushlin  lini  1068°K  it  the  meltin*  point  of  M0O3  . 
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Now  it  is  necessary  to  estimate  the  free  energy  of  vaporizing  liquid  M0O3. 


**°°3(1)  — *'Mo03ff) 


We  can  estimate  that  AF°15j0  -ofor  the  above  vaporization  process  from  the 

fact  that  the  normal  boiling  point  of  MoO.  is  1.  530*  K  according  to  Brewer  39. 

JhiVtS8tiTn^tlvn  l*  Probably  a  conservative  one  because  Berkowitz.  Inghram. 
and  Chupka  have  shown  that,  at  the  boiling  point,  the  species  (MoO,),.  (MoO,K. 
and  (Moo,),.  are  probably  most  important.  They  indicate  that  the  minomir 
is  present  xn  an  amount  less  than  a  few  percent  at  the  temperature  of  their 
measurements  (near  900*  K).  On  the  other  hand.  Brewer  and  Elliott39*  40 
indicate  that  the  monomer  is  probably  the  most  important  specie.  For  the 
present  calculation  we  make  the  simplifying  assumption  that  only  the  monome 
xs  important.  We  therefore  assume  that  Af°  530«O.  Actually,  if  polymers 

of  MoOj  do  exist  at  the  normal  boiling  point,  then  AFf530  >  O  for  the  above 
reaction. 


We  can  now  estimate  the  free  energy  change  for  vaporization  by  assuming 
Trouton  s  constant  to  be  21.  Therefore  AF14?3  =  400  calorie  /mole. 


Now  combining  the  last  two  reactions  (liquid  MoO 
we  find  for  the  overall,  desired  reaction:  3 


formation  and  vaporization) 


^^1,473°K  ”  ♦  63,000  ♦  400  -  63,400  calorie/mole 

1/3  -  AF/RT 

Since  PsiO  *  PMoOj  .  e 

we  find 


PgiO  “  *’24  x  10” 7  atm. 
pMo03  "  •*  *  10  7  «">• 


and 


Ptotml  -  1-65  *  lO”7  atm. 
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Thi.  last  value  of  that  pressure  is  in  good  agreement  with  the  value  of  1.  1  x 
,  .  atm  ”P°rted  by  Inuzuka  38  xhe  above  di,cu„ion  ,how,  he  wa,  aot 

studying  the  vaporization  of  pure  SiOj  in  a  neutral  atmosphere  and  gives 
1  onal  strength  to  the  data  used  in  the  present  report.  Since  Coughlin 
uses  essentially  the  same  heat  of  formation  data  for  SiO,  .  as  we  used  earlier, 
our  recalculation  of  Inuzuka's  results  further  substantiate  that  choice. 

f°  £.ar,W*  hay*  °"ly  considered  one  reaction  in  SiOj  decomposition  as  indicated 

obtain  h  !  tHlS  SeCti°n  ff  the  rep0rt*  Because  fairly  good  agreement  is 
°‘*al“*d  W1*  experimental  data  in  thi.  manner,  it  appear,  probable  that  thi. 
reaction  is  the  most  important  one.  However,  other  additional  competing 
°f  simultaneous  reactions  may  also  be  important.  These  will  now  be  con- 
sidered  in  the  following  sections. 
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HI.  THERMOCHEMISTRY  OF  THE  REACTION  si0(»)  ^^(a) 


The  equilibrium  between  solid  and  gaseous  SiO  is  of  interest  because  of  the 
possibility  that  SK)(fl)  formation  may  limit  the  vapor  pressure  of  S*°(g)  to 
be  expected  from  the  reactions  discussed  in  the  previous  sections.  Data 
concerning  properties  of  solid  are  very  limited.  According  to  Brewer 
and  Edwards^  ' and  Hoch  and  Johnson,  ^ ^  Si°(s)  unstable  relative  to  Si(s) 
and  SiOjt,)  below  1, 450*  K  and  1,573*K,  respectively.  However,  above 
these  temperatures  the  solid  phase  SiO(s)  does  exist  as  shown  by  X-ray  dif¬ 
fraction  results. 


The  limited  amount  of  thermodynamic  data  that  are  available  for  siO^  are 
also  very  uncertain.  Brewer  and  Edwards®  have  interpreted  the  results  of 
Tombs  and  Welch^  obtained  over  a  mixture  of  si  and  SK>2  as  referring  to 
the  reaction  SiO, ,) si^(g)  in  the  temperature  range  of  (1, 573*  -  1,920*  K) 

They  state  that  for  this  temperature  range  AH  *  58,550  calorie /mole  and 
As  -  25.45  e.  u.  ,  although  they  point  out  that  the  entropy  data  may  be  in¬ 
correct.  If  these  data  are  used  as  stated,  for  free  energy  of  reaction  versus 
temperature,  the  approximate  relation  af  -  58,550  -  T*25.45  la  obtained. 

From  this,  one  obtains  p  s  1  atm  (AF-O  )  at  T  =2,  300*  K. 

Additional  vapor  pressure  data  calculated  from  this  relation  are  plotted  in 
figure  6  by  extrapolating  to  temperatures  above  2,  000*  K.  This  extrapolation 
is  not  considered  too  good  in  view  of  the  entropy  uncertainty. 


One  cam  also  utilize  the  above  data  in  another,  more  correct  extrapolation 
above  2,  0Q0*K.  For  this,  the  following  data  were  taken  from  Brewer  and 
Edwards:  8*  19 


C  -  IO.57  ♦  2.98  x  10”  h*  _  2.72  x  105t“2  ; 

P(c$8?n“e 


^298  (SiO  crystalline) 


—  6.4  e.u.  ; 


and 


*298  (SiO  amorphous) 


•  7.3  e.u 
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With  these  data  we  find  that  s|  g(g  .  31.1  *  3.0  ,.u. 

arorphous) 

Using  our  earlier  result  that  S 2  -  66.O4  e.u.  .  we  find  that  for 

•sro^) 

the  present  reaction  .  34 .94  e#u> 

If  we  now  assume  that  the  earlier  relation  af  _  58t55o  _  25.45  T  is  valid  at 
T  -  2,  000  K,  we  can  fit  the  equation  at  2,  000*  K  with  better  entropy  values  to 
extrapolate  to  higher  temperatures. 

By  using  more  realistic  entropy  values,  we  obtain 

AF  -  77,530  -  34.94  T. 


The  above  relation  corresponds  to  a  heat  of  reaction  of  AH  -  77,530 
calorie/mole  in  the  temperature  range  near  2,000*  K.  These  data  are  plotted 

in  figure  1.  The  line  is  dashed  above  2,000*  K,  because  of  the  uncertainty  in 
the  data. 

Using  the  data  of  Brewer39  on  the  heat  of  formation  of  SiO(s)  (  -  104,600 

calorie/ mole )  >  we  can  estimate  the  heat  of  the  reaction,  siO,  r  siO 

as  follows:  (s) 


AH298  ”  AH,  “  AH,  -  —  21,800  —  (404,600)  calorie/mole 

S,°(g)  S,°(s) 


“  +  82,800  calorie/mole 


Also, 


A Ht  -  AH- 


'T  "  ^n298  +  |  ACpdT 

'298 


1 

/ 


AH 


T  “  A H 293  +  (Ht  -  H298)si0  ~  <HT  ~  H298> 


SiO, 


(s) 
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•  *"»•  +  <HT  -  "o’siO^  ~  <H»8  -  H0)s.0(j)  -  («T  -  H298)s.0(>) 


Therefore, 


AH2,000  •  82,800  ♦  8.18  x  2,000  -  6.955  *  298 


-  82,800  ♦  16,360  -  2,070  -  23,090 


AH2,  ,000  "  74,000  calorie/mole. 


This  value  agrees  well  with  the  value  of  77,  530  calorie/mole  given  above, 
and  thus  makes  the  relation  seem  more  reliable.  Accordingly,  the  data  for 

the  relation  AF  -  77,530  —  34*94  T  have  been  taken  as  the  best  and  have  been 
plotted  in  figure  6. 

These  data  lead  to  a  larger  vapor  pressure  of  siO(  )  in  equilibrium  with  siO, 
than  predicted  by  the  first  expression  Af  -  58,500  -  25.45 T  .  However,  the 

data  show  that  the  vapor  pressure  of  SiO^j  ,  based  on  either  expression,  is 
much  greater  than  the  vapor  pressure  of  SiO(g)  in  equilibrium  with  si°2(*) 
under  neutral  conditions.  It  would  ,  therefore,  not  be  expected  that  SiO(#) 
could  be  formed  from  siO_.  .  under  neutral  conditions. 
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iv.  ASSOCIATION  REACTIONS 


Very  little  information  is  available  regarding  the  formation  of  dimers  and 
trimers  of  SiO(  ,  .  The  paper  of  Porter,  Chupka.  and  Inghram  3  gives  one 

measurement  oT  the  si202  pressure  in  equilibrium  with  SiO,  v  .  These  pres- 

sures  are  4.  5  x  10”  atm.  and  1 .  1  x  10“4  atm.  ,  respectively,  at  1, 463*  K. 

From  this  it  follows  that  K  -  ^2° 2  =  3.  72  at  1,463*  K. 

*  2 
PSiO 

From  the  above  equilibrium  constant,  we  can  calculate  that  the  free  energy 

change  for  dimerization  is  Ac  ..  a  o->o  i_  •  /  ,  _ 

_  ®  1.463  -3.820  calorie/ mole  Ot  Si-,0,/  , 

To  get  a  very  rough  estimate  of  She  degree  of  dimerization  at  high2er(fempera 

tures,  it  is  necessary  to  estimate  the  entropy  change  for  the  above  reaction. 

Data  are  already  available  for  the  entropy  of  the  monomer,  siO,  ,  .  For  the 

dimer,  however,  no  molecular  data  are  available,  and  its  entropy  must  be 

estimated  by  approximation  methods. 


The  empirical  equation  of  Kubaschewski  and  Evans *5  can  be  used  to 
estimate  s298  -  64-1  e.t 


for  Si202(g) 


Assuming  that  the  Si,0 


about  the  same  entropy  increment  as  another  four-atom  molecule,  C1F* 

ttlk/,  A  »  — •  A  •  9  J 


has 


whose  entropy  increment  is  s 


-  S 


o  ___  4i  ‘  ' ,  :  '  ^l,5«o  ~  *^298  "  29.J  e-u.  (according  to  the  NBS 

Series  III  tabulation),  we  get  Si,500  -  93.6t3.0e.u-  for  SijO,/ ■  }  ,  where 

the  uncertainty  is  only  estimated.  From  the  tables  of  Edward  JvWe  find  that 

slt500  “  ^3*3  e.u.  ;  and  hence  for  the  reaction  as  given 

Sl°(g) 


AS  -  93.6  -  2(63.5)  -  -33.4  e.u. 


Assuming  that  As,  .qo  -  As, 
A  Ft  —  AH,  -T  x  if 


.  and  using  the  relation 

1 1,463  ”  1  *  aSl,463,  we  find  that  AH,  46,  -  -52,500  calorie/moje  Si, O,,,  . 

From  this  we  find  that  aFt  .  _  52,500  ♦  33-4  T  .  This  relation  is,  of 

course,  very  rough  and  at  an  extrapolated  temperature  of  3,  000*  K  the  un¬ 
certainty  may  be  ±  5,000  calorie/mole  or  more.  The  free  energy  change 
for  the  reverse  reaction  is  shown  in  figure  1. 


At  3,000*  K,  a  free  energy  change  of  +47,700  calorie s /mole  is  obtained  from 
the  above  equation,  corresponding  to  an  equilibrium  constant  of  3.46  x  10““*. 

The  expected  pressure  ofsi202  as  a  function  of  SiO  pressure  is  presented 
in  Table  VIII. 
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TABLE  VIII 


PRESSURES  OF  S»202(f)  in  EQUILIBRIUM  WITH  SiO(g)  AT  3,  000  K 


rSiO  (son) 

Si2°2  (stm) 

0.4 

5. 55  x  10~5 

1 

3. 46  x  10"4 

10 

3.46  x  10“2 

100 

3. 46  atm. 

The  data  in  Table  VIII  shows  that  the  dimerization  effect  is  probably  not  too 
important  for  the  case  where  the  boiling  point  is  above  3,  000*  K  and  the  am¬ 
bient  pressures  are  less  than  50  atms.  However,  the  dimerization  phenome¬ 
non  would  be  more  important  (for  a  given  pressure  of  siO  )  at  lower  tempera¬ 
tures.  Again,  it  should  be  emphasized  that  these  data  do  not  nave  high  ac¬ 
curacy. 
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v.  OXYGEN  DISSOCIATION 


\ric°o™pig.-r  teffrtive 

rangl’o^t  90°0 ‘f’j  toO^K  ha'.T *  °f  ^  dis#ociation  «*  °2  in  the^e^raTure 
been  used  in  the  pr’esent  case  ,t  '  The  data  °£  Stil  and  Sinke^have 

data  from  3,  000*  K  to  3  100*  K  to  oht  necessary  *°  extrapolate  the  free  energy 
temperature  in  °  '  ,  K  *°  obtam  an  «9“ilibrium  constant  at  the  higher 

da"aP7  xL  Halils  of  fhe^o  Kh  ?***  Wi*h  th*  earli"  NBS  Series  m  4 f 

values  of  the  equilibrium  constants  obtained  are  shown  in  Table  IX. 


EQUILIBRIUM  constants  for  THE  REACTION  X°2(g) 


>■*) 


Temperature 

CK) 

2,900 
3.  000 
3,  100 


.  po 


0.  079 
0.  112 
0.  157 


From  the  equilibrium  constants  shown  in  TaHl*.  tv  ,  ,1  ,  . 

ratio  have  been  obtained  for  various  selects  *  V  UCS  °f  the  pOy/Po 

results  of  these  calculations  p^d** £1*"-*”“*  ^ 

rroo^KiuTheis  alsn  :bout  °- 01  a-  0°z?iz  in 

3.  100  K,  there  will  be  an  atomic  oxygen  pressure  of  .  ’  ‘ 


2 

itm 


Sr  ir 

get  an  approximate  linear  relation  for  the  free  *  dl8soclat,on  reaction.  To 
at  3,000‘K,  the  data  of  Stull  and  Sinke25was  used  'fIo  °  .th'  reaction 

ThT.Tedlatir°Py  data’/he  f°ll0Win*  relation  is  obtained  A  ,-«.£*  °\lnr' 
rbxs  relation  gives  a  free  enercv  chancre  of  ?q  i  ?o  i  .  ,  61,348  -  16. J 1  T. 

whereas  Stull  and  Sinke  •  s^tabiultion  1,  29  091  ,  **  00°*  K- 

free  energy  relation  and  the  '  alorie/mole.  The  above 

in  figures  [  and  2  corresponding  equilibrium  constant.  are  plotted 
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The  data  on  the  dissociation  of  oxygen  are  well  known  and  provide  a  test  of 
the  reliability  of  the  SiC>2  decomposition  data  reported  by  Porter  et  al,  inas¬ 
much  as  their  mass  spectrometric  measurements  included  both  Oi  and  O 
At  1, 900*  K  they  give  2 


PO  -  2.2  m  10  “7  atm.  and  p0^  -  8.2  *  10  atm. 
for  the  reaction  .  This  corresponds  to 

*1,900  -  7.7  *  10-4  -  pQ 

Po2^ 


Hence  AF.  jqq  27,065  calorie  /mote 
Stull  and  Sinke**'*  give 


F or  the  same  dissociation  reaction, 


Temperature,  *  K  LogipK 


AF°,  calorie/mole 


’.900 
2,  000 
2.  100 


-3. 529 
-3. 178 
-2.  859 


30, 678 
29, 091 
27, 477 


Thu.,  we  lee  that  Porter',  free  energy  value  of  27,065  calorie/mole,  based 
in900**KP°rted  rr***Ure  meaaurement»,  is  low  by  3,613  calorie/mole  at 

“  .h*d  *tternpted  to  correct  Porter's  vapor  pressures  by  using  the  more 
recent  gold  vapor  pressure  data,  then  the  pressures  would  have  been  increased 
by  a  factor  of  3.  95.  This  would  increase  the  equilibrium  constant  at  1, 900*  K 

to  a  vaiue  "f  lo810K  =  -2-  815,  thus  making  the  standard  free  energy  change 
1,900  *  24,400  calorie/mole  .  ** 

However  this  value  is  in  error  by  6,  300  calorie,  at  1, 900*  K,  whereas  the 

milU„  V™Ue  r,eported  ^  port«r,  etal,  i.  i„  error  by  only  3,600  calories  / 
mole.  Therefore,  it  seem,  that  there  may  have  been  either  compensating 
error,  in  Porter's  reported  data,  in  their  choice  of  cross  sections,  or  the 

uronef^s  moPOr,  preVUr'  data  may  be  in  m  any  case,  at  this  time 

it  appear,  more  logical  to  accept  Porter's  tabulated  pressure  data  as  being 

correct  because  it  yield,  free  energy  values  in  better  agreement  with  the 
known  values  on  the  02  dissociation  reaction. 
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VI.  THERMOCHEMISTRY  OF  THE  REACTION  ^(g)-*  ^(g)  *  ‘^(g) 


To  estimate  the  importance  of  the  siO  dissociation  reaction  to  yield  sif  \  we 
can  use  the  following  free  energy  and  entropy  of  formation  data  at  3,  00u*  K: 


For  Si  (g),  AFj(000  "  10,618  calorie/mole  (Ref.  25,  Stull  a.  d  Si  nice) 

using  Si^  j j  as  reference  state, 


ForOjCg),  AF}iooo  “  ®  (Since  Oj  is  reference  state). 


For  the  formation  of  SiO,  }  (i<e.t  for  the  reaction 
Si(l)  +  /i02(g)  -*■  Si 


A F  jtooO  “  ^H3,000  -T  ASjtooo  "  -38,954  -  3,000(ASj  Q001  calorie/mole 


As3,000  “  SsiO  ~‘^S02(g)  “  ^‘(l)  “  10-17  e‘U’ 


and 


AF3  000  “  69,464  calorie/mole  (for  SiO(gj  formation) 


where  the  entropies  of  Si^  and  are  from  Stull  and  Sinke^and  that  of 
from  our  calculations.  Assuming  that  AH  and  AS  do  not  depend 


SiO, 


(g) 

upon  temperature,  we  find  that  AF  »  -  38,954  -  10.17  T  for  this  reaction  in 
the  temperature  range  of  2,  000°  K  to  3,  000*  K.  From  these  data,  we  find 
that,  for  the  first  reaction  given  above,  (i.e.  si0(g)  “ *  s*(g)  *  ^°2(g)  ) 


AF3000  ”  +  80,082  calorie/mole 
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and 


As3,000  "  *  s02(k)  *  Ssi(g)  “  ^S»°(g)  “  16,31  e,u* 

Hence.  AF  -  129,000  -  16.31  T  near  T  -  3.000°K  , 

Now  we  can  calculate  the  equilibrium  constant  at  3,  000*  K, 


AF3  ooo  *  “  RT  In  K  -  90,082 

„  j'  PsiPO, 

or  K  -  1.45  *  10“ 6  - - 2_ 

pSiO 

If  we  consider  the  case  where  PsiO  “  °*4  son.  and  Po2  “  °*2  atm.  t  (i.©., 
the  approximate  vapor  pressures  under  neutral  conditions  at  3,  000*  K  from 
figure  3)  we  find  that 

pSi  -  1.30  x  10  “6  atm. 


It  can  be  seen  that  a  negligible  amount  of  monomeric  S*(g)  vapor  exists  under 
these  conditions.  Dimers  and  larger  aggregates  of  Si  would  also  be  negligi¬ 
ble  under  these  conditions  since  the  equilibrium  studies  of  Honig  24  show  the 
monomer  to  be  the  most  important  species,  even  when  the  condensed  phase  of 
silicon  is  present.  The  calculated  partial  pressure  of  silicon  is  far  removed 
from  such  saturation  conditions  inasmuch  as  Stull  and  Sink^^  report  a  boiling 
point  of  Silicon  of  Z,  950*  K  and  Honig  ^  gives  3,  060*  K. 
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VII.  THERMOCHEMISTRY  OF  THE  REACTION  SiC^p  — » Si02(g) 

A.  GENERAL 

For  estimation  of  thermochemical  quantities  involved  in  vaporization  of  liquid 
Si02  to  gaseous  Si02  ,  the  most  pertinent  available  data  are  those  of  Porter, 
Chupka,  and  Inghram.  ^  They  have  made  a  mass  spectrometric  study  of  the 
decomposition  of  cristobalite  at  1,  800*  K  and  1,  900*  K. 

S‘®2(cri«tob*lice) 


They  found  that  =  ±  ®  Kcalorie /mole,  from  slopes  of  their  log 

i/T  data.  However,  in  a  personal  communication  Porter  has  stated 
that  more  accurate  heat  of  reaction  data  (and  hence  free  energy  data)  could 
be  obtained  for  the  above  reaction  with  a  suitable  estimate  of  the  entropy 
change  in  combination  with  their  published  pressures.  Brewer  and  Searcy ^ 
have  also  discussed  the  value  of  this  type  of  calculation. 

B.  ESTIMATE  OF  ENTROPY  FOR  Si02(g) 

Since  there  is  little  published  information  about  Si02^g)  ,  it  has  been  neces¬ 

sary  to  make  some  assumptions  about  the  configuration  of  the  molecule. 
Knowledge  about  the  c02  molecule  is  used  as  the  basis  for  estimating  proper¬ 
ties  of  the  SiO  2  molecule. 


Since  the  c02  molecule  has  a  linear  structure,  it  has  been  assumed  that  Si02 
is  also  linear.  In  estimating  the  interatomic  distances  the  following  data 
were  considered: 


Molecule 


SiO 

CO 

co2  -  o> 


f*  (A? 
1.  510 

1.  1284 

1.  1615 


Ref. 

U  u  32 

Herzberg 

43 

He  rzberg 

u  .  43 

Herzberg 


From  these  data  we  estimate  the  interatomic  distance  (si-O  )  of  linear  SiO 
to  be  1.  54A. 


The  estimation  of  the  vibrational  frequencies  for  the  linear  si02  molecule 
has  also  been  done  by  analogy  withcOj  .  In  Table  X  are  shown  some 
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vibrational  frequency  data  for  Ct>2  ,  cs  .  CO  .  and  cs  molecules.  It 
is  seen  that  for  the  diatomic  molecules  the  ratio  of  the  vibrational  frequencies. 

res/ -co  =  °-  59  which  agrees -qualitatively  (maximum  error  =17  percent) 
withthe  corresponding  ratios  for  the  triatomic  molecules  c02  and  cs,  .  The 
corresponding  ratio  of  vibrational  frequencies  for  the  diatomic  SiO  ana  CO 
was  obtained  from  Herzberg's  book  3*which  gives  vSiO/vcO  =  0.572.  This 
ratio  was  then  used  to  calculate  the  vibrational  frequencies  of  si02  from 
the  corresponding  frequencies  for  C02  •  This  calculation  is  approximate 
and  the  calculated  vibrational  frequencies  might  be  in  error  by  as  much  as 
15  to  20  percent. 


TABLE  X 

VIBRATIONAL  FREQUENCY  DATA  ESTIMATED  FOR  SK>2(5) 


Mode 

co2 

cs2 

Ratio  vCS^/vCO? 

Ref. 

Stretching. ul 

1 .  342  cm”  * 

656.  5 

0.  49 

Herzberg^^ 

Bending.  v2 
(degenerate) 

667 

396 

0.  59 

Herzberg43 

Stretchi  .^,, 

2,  349 

1. 523 

0.65 

Herzberg43 

CO 

CS 

Ratio  vCS/vCO 

*1 

2,  170 

1, 285 

0.  59 

Herzberg^ 

CO 

SiO 

Ratio  vSiO/vCO 

vl 

2,  170 

1. 242 

0.  572 

Herzberg3^ 

co2 

Si02 

Ratio  uSi02/vC02 

(calculated) 

(assumed) 

v  1 

1.342  cm"1 

767 

0.  572 

vl 

667 

382 

0.  572 

*3 

2.  349 

1.  340 

0.  572 
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For  the  moment  of  inertia  of  the  si(>2  molecule,  we  used  the  relation  given 
by  Herzberg  for  the  moment  of  inertia  of  a  symmetrical  molecule. 

*B  “  2MorSi  -  O 

where  Mc  is  the  mass  of  the  oxygen  atom  and  rsi  _0  is  the  interatomic  dis¬ 
tance 

M0  *  16  z  1 .66  x  10  ~  24  gm 

fSi  -  O  "  M4A°  -  1.54  x  10“8  Cm, 

Hence,  we  estimate  for  siC>2  that 
Ig  “  1.26x10“  ^8  gm  — cm2 


For  the  temperature  ranges  of  possible  interest  in  this  calculation  ( T  <5,  000*  K), 
electronic  contributions  are  assumed  to  be  negligible.  It  is  assumed  that 
only  a  single  electronic  ground  state  is  occupied  because  the  si02  molecule 
has  an  even  number  of  electrons.  Dr.  Beckett,  of  the  National  Bureau  of 
Standards  in  a  private  communication  has  stated  that  very  roughly  speaking 
the  n  higher  state  would  probably  be  at  least  5  e.  v.  above  the  ground 
state.  Since  5  e.  v.  corresponds  to  a  value  of  RT  =  115,300  calorie,  or 
T  =  58,  000*  K,  it  can  be  seen  that  contributions  to  the  partition  function  will 
be  negligible  at  the  temperatures  under  consideration  here  and  essentially 
only  the  ground  electronic  state  will  be  occupied. 

The  calculation  of  the  entropy  for  si°2(g)  was  ma<*e  with  the  statistical 
mechanical  formula  given  by  Kelley14  for  linear  polyatomic  molecules. 


S  m  3/2  R  ln^M  +  7/2  R  lneT  —  R  In  P  -f  R  In  I  —  R  lna  —  +  2S  (Einstein)  4-  ♦  17J.385 


The  present  calculation  is  for  atmosphere  pressure  and  so  the  third  term 
above  does  not  contribute  to  S.  The  symmetry  number  <7  has  the  value  two  for 
the  Si02  molecule.  Evaluation  of  the  2  S(Ein»tein)  contribution  must  be 
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made  for  the  four  vibrational  degrees  of  freedom  for  a  linear  triatomic 
molecule. 

Terms  for  ^  and  ^  contribute  once  each  to  this  sum  whereas  those  for 
u2  (the  bending  vibration)  contribute  twice.  The  tables  of  Taylor  and 
Glasstone44  have  facilitated  evaluation  of  the  vibrational  terms.  No 
attempt  has  been  made  to  incorporate  refinements  such  as  Fermi  resonance 
which  apply  to  the  C02  molecule. 

Enthalpy  calculations  were  made  with  the  relation 

(H  -  Hp)  -  7/2  R  -  6.954 

T 

Translation  ♦  rotation  4  el 


The  terms  due  to  vibration  can  again  be  added  from  Taylor  and  Glasstone 's  44 
tables. 

In  Table  XI  are  listed  the  resulting  heat  capacity,  enthalpy,  and  entropy  data 
calculated  for  the  linear  Si02(g)  molecule. 


TABLE  XI 

HEAT  CAPACITY,  ENTHALPY,  AND  ENTROPY  DATA 


FOR  THE  ASSUMED 


Temper¬ 
ature  (•  K) 

cp 

(calorie  / 
mole  *  K) 

(H-Hq) 

T 

(calorie  / 
mole*  K) 

298 

10.  80 

8.  53 

500 

12.  55 

9.  83 

1. 000 

14.  10 

11.  66 

1,  500 

14.  50 

12.  57 

1, 800 

12.  90 

1. 900 

12.  99 

2.  000 

14.  67 

13.  06 

2,  500 

14.  73 

13.  40 

3.  000 

14.  79 

13.  62 

LINEAR  si°2(g)  MOLECULE 


H-H0 
(calorie  / 
mole) 

H-H298 
(calorie  / 
mole) 

S 

(e.  u.  ) 

2,  542 

0 

54.  65 

4,  915 

2,  373 

11, 660 

9,  118 

18,  855 

16, 313 

23,  220 

20,678 

78.  57 

24,  681 

22,  139 

79.  37 

26,  120 

23,  578 

80.  06 

33,  500 

30,  958 

83.  37 

40,  860 

38,  318 

85.  96 
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The  heat  capacity  data  obtained  in  Table  XI  are  plotted  in  figure  8.  For  com- 
panson  purpose*  the  heat  capacity  of  CO 


2  (*) 


has  been  included.  It  is  seen 
exceeds  the  estimated  heat 


that  above  2 ,  500*  K  the  heat  capacity  of  cO 

capacity  for  SiC>2(g)  .  This  probably  occurs  because  such  complicating 
factors  in  c02  as  the  Fermi  resonance  (or  accidental  degeneracy)  and  other 
refinements  have  not  been  considered  in  the  present  calculation. 


From  Porter  et_al's  paper  3  ,  we  obtain  for  the  reaction  siO.,,  .  x  siO 
the  vapor  pressure  data  shown  in  Table  XII.  The  free  energies"  sLwn  we& 
calculated  from  their  vapor  pressure  data. 


TABLE  XII 

DATA  FOR  Si02(cfi,t)— *si02(g) 

1 ,  800*  K  1 , 900*  K 

Psi02  (atm.)  1.6  x  1 0“8  l.OxlO"7 

AF  —  RT  In  psio2  +64,  214  +60,860 

(calorie /mole) 


sSi02(g)  (from  Table  XI)  e.  u. 

Ssi°2(cri«t) 

(from  Kelley1*9)  e.  u. 

A  S  e.  u. 

T  A  S  calorie  /mole 
AHt  -  AF  ♦  TAS 
calorie /mole 


78.  57 
38.  15 

40.  42 
72, 756 
1 36, 970 


79.  37 

39.  12 

40.  25 
76, 475 

137, 335 
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1°  *<»8  *  reaction  data  at  1 . 800 ‘and  1. 900*  X  to  heat  of  reaction 

’  .K*  we  can  “•*  **»e  Jalue»  from  Table  XI  for  the  heat  content  of 

way  one  obtai  K«“«»I3fos.that  of  cri.tobailite  at  1,  800-  K.  In  this 


*”298  "  *”t  -  (Ht  -  H298) 


Si  O 


2(g) 


(Ht  -  H298) 


SiO 


2(crist) 


*^298  “  140,822  calorie/mole 


Also,  from  data  at  1, 900*  K,  we  obtain 
*”298  “  141.516  . 


The  average  of  the.e  two  result,  is  AH29g  =  141.  169  calorie  /mole.  This 

H  t;  !Toooat  tha"  the  °riginal  value  reported  by  Porter  et  al  of 

136.000  ±  8,000  calone/mole  but  falls  within  their  range  of  uncertainly 


*”2.000 

follows : 


for  the  reaction  SiO 


2(crisc)‘ 


SiO 


2(g) 


,  is  calculated  as 


AH 


2,000  -  *”298  ♦  <”T  ”  H298)  c- 


SiO 


-  (Ht  - 


2(g) 


T  “  H298>  • 


criat.. 


*”2,000  "  136,627  calorie/mole 


By  taking  into  account  the  heat  of  fusion  of  cristobalite  (2,  500  calorie /mole) 
we  can  calculate  the  heat  of  vaporization  of  liquid  silica 


SiO 


2(1)  ^  ^*®2(g)  *”2,000  "  134,127  calorie/mole 


from  the  above. 
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For  the  heat  of  this  reaction  at  higher  temperatures  we  have 


AH-r  -  AH 


2,000 


+  (Ht  -  H 


2,000* 


SiO 


-  (H-r-  _  H 


2(g) 


2,000*, 


SiO 


2(1) 


The  additional  necessary  data  for  calculations  of  AHt  are  in  Table  I  and 
a  e  !•  Results  of  such  a  calculation  are  summarized  in  Table  XIII. 


TABLE  XIII 

HEAT  OF  VAPORIZATION  OF  LIQUID  si°2 

Si°2(  1)  *’Si°2(g) 

2, 000*  K  2,  500* K  3,  000*  K 


(Ht  -  H2  000) 

Si02(i) 

0 

10,  830 

21, 660 

(calorie  /mole) 

(Ht  -  »2,000*  50 

Si°2(g) 

(calorie  /mole) 

0 

7,  380 

14,  740 

AH 

(  calorie /mole) 

134, 127 

130, 677 

127, 207 

Likewise,  to  calculate  the  entropy  change  for  the  last  reaction. 
Table  VI  and  Table  XI  can  be  used  as  shown  in  Table  XIV. 


the  data  of 


The  free  energy  data  given  in  Table  XIV,  have  been  plotted  in  figure  1.  where 

it  can  be  seen  .t  AF  «  O  (PsiO-%  ”  1  *tm**af  rr_  * 

aiw2  at  T-  3,  540*  K.  These  data  show 

that  at  3,  000*  K,  Af  -  19.416  calorie /mole  and 


-  19,416 

Pcjo  "  e  — - - 

blU2  Rx3,000 


O.O39  utm. 
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TABLE  XIV 


ENTROPY  CHANGE  FOR  VAPORIZATION  OF  LIQUID  Si02 


StO 


^(l) 

AS 

tas 

AH 

AF  -  AH  -  TAS 


(e.  u.  ) 

(e.  u.  ) 

(e.  u.  ) 

(calorie  /mole) 
(calorie  /mole) 
(calorie  /mole) 


2(1) 


SiO 


2(f) 


2, OOP*  K 
80.  06 


2, 500* K 


83.  37 


3, 000*  K 
85.  96 


41.  07  ±  0.6 


46*  13  ±  0.  8  50.03  ±1.0 


38.99 


77, 980 
134,  127 
56,  147 


37.24 
93,100 
130,677 
37, 577 


35.93 

107,791 

127,207 

19,416 


SiOjfg)  pressures  as  a  function  of  1/T  are  shown  in  figures  2  and  11.  In 

this  case  the  equilibrium  constant  is  equal  to  the  partial  pressure  of  Si02(s) 

in  atmospheres.  If  si02  vaporizes  according  to  the  above  reaction  at 

3,  000*  F,  then  the  total  amount  of  heat  absorbed,  including  the  sensible  heat 
of  solid  and  liquid  Si02  from  296*  to  3,  000*  K,  is 


"  Hj98>a°2(i>  ‘  AH<Y*^mon) 

-  <H2,000  -  H»s)Si02(gU..)  +  <H3»®°°  "  H2,00°>sio2(1*  A  Hr^>oril«ion 
"  177,877  c*lorie/mole  (5,336  Btu/lb) . 
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The  sensible  heat  (the  first  two  terms  in  the  second  equation)  contributes 
50,670  calorie/mole  (1,520  Btu/lb.  ),  and  the  vaporization  contributes  127,207 
calorie  /mole  (3,  816  Btu/lb).  It  should  be  noted  that  simple  vaporization  of 
Si02(s)  at  3»  000#  K  gives  a  total  heat  absorption  of  only  5,  336  Btu/lb  as  against 
the  valueof  6,650  Btu/lb  found  for  the  case  of  SiO(g)  and  02fg)  formation. 
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VIII.  SiOz  IN  REDUCING  ATMOSPHERES 


In  a  reducing  atmosphere  decomposition  of  Si02,  to  produce  SiO^j,  occurs 

at  a  much  lower  temperature  than  under  neutral  or  oxidizing  conditions. 

Three  typical  reducing  agents,  for  which  data  are  available,  include  elemen¬ 
tal  silicon,  carbon,  and  hydrogen.  Since  the  reducing  reactions  occur  at 
lower  temperatures  they  are  of  possible  interest  for  at  least  two  important 
reasons.  First,  a  re-entry  vehicle  might  be  constructed  out  of  a  mixture  of 
Si02  and  a  reducing  agent.  This  might  permit  decomposition  of  Si02  to  occur 
at  a  lower  temperature  where  its  viscosity  is  still  high.  Second,  knowledge 
of  the  effect  of  reducing  reactions  is  important  because  many  of  the  current 
testing  procedures  in  arcs  are  being  carried  out  under  reducing  conditions. 
That  is.  arc  testing  done  with  carbon  electrodes  which  yield  carbon  atoms  in 
the  plasma  may  really  correspond  to  a  reducing  atmosphere  when  carbon  is 
present  in  excess  of  oxygen  or  if  the  other  ingredients  of  the  plasma  are 
present  in  stoichiometric  amounts  of  reducing  and  oxidizing  agents.  Also, 
if  a  test  specimen  in  the  arc  is  subject  to  stoichiometric  amounts  of  H  atoms 
and  OH  atoms  (reducing  agent  and  oxidizing  agent  respectively),  the  higher 
mobility  of  the  H  atoms  may  effectively  make  the  specimen  "see"  a  reducing 
atmosphere. 

The  reaction  between  stoichiometric  amounts  of  silicon  and  silica  leads  to 
the  reaction: 

Si  +  SiO^-^2  SiO(g) 


Although  Si  +•  Si02  in  a  closed  system  will  decompose  at  a  lower  temperature 

than  either  material  separately,  it  cannot  be  concluded  that  the  reaction  is 
independent  of  oxygen  pressure.  This  arises  from  the  fact  that  the  simultaneous 
reaction  Si  +•  02  must  be  considered.  In  the  dynamic  case  involving  decomposi¬ 
tion  of  the  re-entry  vehicle  surface,  it  is  possible  and.  very  likely,  that  the 
reaction  Si  4*  Si02  will  be  the  most  important.  However,  if  one  is  discussing 
an  equilibrium  case  where  an  unlimited  amount  of  oxygen  is  available,  then 
the  process  reverts  to  a  decomposition  of  Si02  in  an  oxidizing  atmosphere. 

Some  typical  vapor  pressure  data  for  the  Si-Si02  reaction  taken  from 
Humphrey,  Todd,  Coughlin,  and  King18  are  plotted  in  the  lower  curve  of  figure  9. 

For  the  reducing  reaction  between  carbon  and  silica,  one  of  the  possible 
reactions  is  that  for  the  formation  of  SiO^g)  and  CO^gj. 
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si°2(a)  +  ♦  co(g) 

» 

The  data  for  this  reaction  were  taken  from  Humphrey's  tables*®and  plotted  in 
figure  9.  The  upper  curve  gives  the  CO  pressure  while  the  lower  curve  is 
again  assumed  to  represent  the  actual  SiO  pressure  in  the  system  because  of 
disproportionation  to  Si|s|  and  SiO^^p  Reference  to  these  data  would  indicate 

that  a  pressure  of  one  atm  (  p^-q  -  0.89  »tm,  p^q  “  0*11  *tm  )  reached  at 

a  temperature  of  1, 845*  K.  The  assumption  that  the  SiO  pressure  is  fixed  by 
the  Si-SiO^  equilibrium  appears  to  be  quite  valid  since  if  it  had  been  ignored 
higher  SiO  pressures  would  have  been  obtained  by  calculation.  Of  course, 
the  SiO^Bj  =  SiO^gj  equilibrium  instead  of  the  Si-Si02  equilibrium  may  fix 

the  SiO^gj  pressure  in  the  present  case.  Actually  these  two  equilibria  are 

probably  quite  similar  (i.  e.  ,  give  the  same  SiO^gj  pressure  at  temperatures 

near  1.  700*  to  1.800*  K)  if  we  assume,  as  discussed  in  section  II.  that  SiO(8) 

becomes  stable  near  1,450*  to  1,  573*K. 

Another  reducing  reaction  of  interest  is  that  of  SiO^  with  H^.  The  applicable 
experimental  data  of  Grube  and  Speidel21,  as  given  by  Edwards1^  have  been 
plotted  in  figure  10.  Here,  it  can  be  seen  that  their  tabulated  vapor  pressure 
data  for  SiO^g)  are  in  agreement  with  those  determined  by  Tombs  and  Welch^. 

Figure  10  also  shows  expected  vapor  pressures  of  SiO^gj  corresponding  to  the 

vaporization  of  Si0^sj  or  Si^,j  +  ^^2(s)*  ^  should  be  noted  that  in  an  H^ 

atmosphere  the  pressure  of  SiO(g)  is  less  than  either  of  these  equilibria  would 

predict  and  hence  SiO  gas  will  be  formed  rather  than  a  solid  phase.  These 
results  show  that  decomposition  of  SiO^  under  these  conditions  occurs  at  much 
lower  temperatures  than  under  oxidizing  or  neutral  conditions. 
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IX.  SUMMARY  OF  CALCULATED  VAPORIZATION  DATA  FOR  SiO mi 

AND  COMPARISON  WITH  PUBLISHED  VALUES 


In  the  preceding  sections  various  equilibria,  involving  SiO£,  have  been  con¬ 
sidered.  The  data  show  that  the  first  reaction,  Si°2(l)^®*°(g)  +  ^  °*(g) 

is  most  important  for  neutral  conditions.  However,  the  reaction  Si°2(ir~>s*02(g)» 
makes  a  significant  contribution  to  the  vapor  pressure.  It  becomes  even  more 
important  under  oxidizing  conditions. 


For  each  of  the  reactions  considered,  a  standard  free  energy  of  reaction  has 
been  obtained.  Simple  approximate  formulas  have  been  developed  to  express 
these  free  energies  as  functions  of  temperature.  The  resulting  formulas  are 
summarized  in  Table  XV.  From  these  approximate  relations  or  the  some¬ 
what  more  accurate  data  obtained  earlier  in  this  paper  for  a  few  of  the  reac¬ 
tions,  equilibrium  constants  have  been  calculated  by  means  of  the  relation 


K 


“  ^f°/rt 

«■  e 


The  results  are  plotted  in  figure  2. 


Although  vaporization  phenomena  were  discussed  in  terms  of  the  single  reac¬ 
tion,  SK>2— H  °2(g)  •  if  is  now  possible  to  consider  an  addi¬ 

tional  reaction,  the  dissociation  of  oxygen.  The  effect  of  the  latter  can  be 
evaluated  both  for  the  case  of  neutral  conditions  and  for  the  case  of  oxidizing 
conditions. 


A.  NEUTRAL  CONDITIONS 

For  the  case  of  neutral  conditions  we  start  with  the  reactions 


Si°2(l)  <  Si°(g)  1/5  ° 2(g)  K1  . 

*02(,)?°<,>  K2  , 


and  let 


PsiO 

PO, 


partial  preaaure  of  SiO, 
partial  pressure  of  O2, 


(8) 

(9) 
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FREE  ENERGY  RELATIONS  FOR  VA 


Pq  m  partial  pressure  of  O. 

Under  our  assumed  conditions ,  we  have  the  additional  relations 

(10) 

(11) 


p02  *  Kp0  "  ^PSiO» 
K1  "  PSiO  *  pOj  • 


(12) 


The  above  three  equations  are  sufficient  to  determine  the  three  unknowns  in 
terms  of  the  two  equilibrium  constants,  P^q  »  po2  »  sndP0. 

Rearranging  equation  (12)  to 

pO-K2pO*  (13) 

and  substituting  into  equation  (10),  we  obtain, 

2  p02  ♦  k2P0^  "  PSiO  (14) 

Finally,  we  substitute  (14)  into  (11)  and  obtain  the  cubic  equation, 

2  P<^2  +  K2  POj  ~  K1  "  °  05) 

for  Po2  in  terms  of  Ki  and  K2  . 

For  the  case  where  k2  is  small,  i.  e.  ,  where  the  dissociation  of  is 
negligible,  we  can  solve  equation  (15)  readily  to  obtain 

P0j  •»  (K,/2)2/J  (K,<<K2P02)  (16) 
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For  the  case  where  we  cannot  neglect  dissociation  of  0^»  equation  (15)  must 
be  solved. 

After  solving  equation  (15)  for  Po£  we  can  obtain  Pq  from  equation  (13)  and 
Psio  from  equation  (10). 

From  selected  values  of  the  equilibrium  constants  Kj  and  at  2,  000*  ; 

2,  500*  ,  3,  000*  ,  3,  500*  ,  and  3,  7  50*  K,  it  has  been  possible  to  calculate 
corresponding  values  of  the  partial  pressures  of  SiO,  O^,  and  O  as  described 
above.  The  resulting  calculated  partial  pressures,  based  on  neutral  condi¬ 
tions  and  including  dissociation  of  oxygen,  have  been  plotted  in  figure  11. 

Referring  to  the  data  in  figure  11,  one  can  see  that  a  total  pressure  of  one 
atmosphere  is  obtained  at  3,  070*  K.  This  corresponds  to  the  normal  boiling 
point  of  SiC>2*  By  reading  partial  pressures  of  the  various  components  from 
the  graph  at  3,  070*  K,  one  obtains  the  following: 


PsiQ  -  0.62  acm 

Pn  •  0.26  atm 
u2 

Pq  —  0.074  atm 
pSi02“  °*P58atm 

"  0.0001  atm 
Ptotal  ”  1*012 


(  +  0.6)  ot  —  0.2)  , 

C  +  0.25  or  —  0.08 )  , 

(  +  0.07  or  -  0.002)  , 

(  +  0.12  or  -  0.002)  , 

(  +  0.0003  or  -  0.00003)  , 
(  +  1.04  atm  or  —  0.28  atm) 


The  uncertainties  in  the  partial  pressures  are  indicated  in  parentheses.  For 
example,  PgiO  may  possibly  have  values  from  0.  42  to  1.  22  at  3,  070*  K. 

It  is  seen  that  within  the  accuracy  of  these  data,  a  total  pressure  of  one  atmo 
sphere  is  obtained  at  3,  070*  ±  75*  K  for  a  neutral  equilibrium  mixture.  The 
accuracy  of  this  boiling  point  can  be  determined  from  the  fact  that  the  free 
energy  function  of  importance  here  is 

AF°  -  173,000  -  54.2  T 


where 
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AF°  is  accurate  to  ±  4*000  calories /mole  near  3,000*  K. 

Using  this  uncertainty  in  AF°  we  obtain  =  54. 2,  from  which 

dT 

it  follows  that  AT  =  ±  75*  K. 

A  more  complete  analysis  of  the  expected  equilibrium  gaseous  species  at 
3,  000*  K  is  as  follows: 


Species 

Pressure  (atm) 

SiO 

0.43 

Si02 

0.  04 

l  si2°2 

5  x  10‘5 

1 

(SiO)x 

* 

unknown 

!  °2 

0.  18 

O 

0.  05 

Si 

1.  3  x  10“6 

Ptotal 

0.  70  atm 

The  results  summarized  in  both  of  the  above  tables  show  that  SiO  and  O2  are 
the  predominant  species.  In  the  neutral,  equilibrium  mixture  Si  and  SigC>2 
are  present  in  minor  amounts.  However,  significant  amounts  of  Si02  and  O 
are  also  present. 

As  we  have  just  seen,  a  normal  boiling  point  of  3,  070*  K  is  predicted  for  Si02* 
If  all  the  other  reactions  had  been  neglected  and  only  the  decomposition  reac¬ 
tion  siOj— ►  SiO(g)  +  KC>2(g)  had  been  considered,  then  a  boiling  point 

of  3,085*  K  would  have  been  predicted,  as  discussed  earlier.  Thus,  the  other 
reactions  have  only  a  small  effect  on  the  predicted  boiling  point.  The  possi¬ 
bility  of  reactions,  other  than  those  already  discussed,  should  not  be  com¬ 
pletely  overlooked,  but  they  are  not  likely  to  be  important.  In  the  absence 
of  further  information  and  in  view  of  the  fact  that  the  observed  reactions  pre¬ 
dict  vaporization  phenomena  in  agreement  with  experiments  at  lower  tempera¬ 
tures,  it  does  seem  likely  that  the  present  calculations  should  give  a  fair 
figure  for  the  boiling  point. 
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Table  XVI  contains  values  of  boiling  points  given  in  the  literature  which  can 
be  compared  with  ours. 


TABLE  XVI 

TABULATED  BOILING  POINTS  FOR  SiOz 


T  ernperature 


Original  Reference 


Date 


Quoted  bi 


2, 230*  C 


Ruff  and  Schmidt^ 


1921  (Lange  -Handbook  of 
Chemistry) 

(Handbook  of  Chemistry 
and  Physics46) 

(Int.  Critical  Tables47) 


2,  590*  C 

5 

Ruff  and  Konschak 

1926 

(Handbook  of  Chemistry 
and  Physics46) 

(Sidgwick48.  Eitel4**) 

2,  800*  C 

Present  work 

2,  950*  C 

O.  Ruff6 

1935 

(A.  Silverman50) 
(Campbell 1 1 ) 

3,  500*  C 

W.  R.  Mott51 

1918 

In  reviewing  the  above  data  it  was  found  that  Brewer3**  has  rightfully 
criticized  the  boiling  point  data  of  Ruff  and  Schmidt4  because  the  measure¬ 


ments  were  carried  out  in  a  carbon  vapor  (reducing)  atmosphere.  The  original 
paper  by  Ruff  and  Konschak  5  also  shows  uncertainty  of  the  same  type.  Al¬ 
though  they  used  an  iridium  container,  it  was  inside  a  graphite  vessel.  Under 
such  conditions  they  observed  rapid  attack  of  the  iridium  container  and  were 


able  to  obtain  only  a  single  measurement  at  2,  060*  C.  They  used  this  measure¬ 
ment  to  obtain  their  extrapolated  boiling  point.  In  view  of  the  fact  that  their 


container  was  subject  to  such  rapid  attack,  it  appears  possible  that  a  reducing 
atmosphere  (of  iridium  or  carbon)  was  actually  present.  This  would  lead  to 
observation  of  a  higher  vapor  pressure  and,  hence,  a  lower  boiling  point  than 
the  true  one.  Other  container  materials  tried  by  Ruff  and  Konschak  5  without 


success  included  tantalum,  tungsten,  molybdenum,  and  zirconium  oxide. 

Ruff  and  Konschak5  further  state  that  their  calculated  boiling  point  of  2,  590#  C 
is  only  a  minimum  value. 
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The  boiling  point  of  2,  950*  C  given  by  Rufi^  in  1935  was  simply  a  restate¬ 
ment  of  his  earlier  work  of  1926.  In  1935  he  stated  that  the  original  2,  590*  C 
boiling  point  was  a  misprint  and  that  the  measured  11.0  mm  pressure  at 
2,  060*  C  should  extrapolate  to  a  boiling  point  of  2,  950*  C.  However,  in  re¬ 
pealing  his  extrapolation  we  find  that  the  2,  590*  C  boiling  point  appears  more 
reasonable.  In  figure  12,  for  example,  the  extrapolation  through  the  ordinate 

11.  1  at—  =0,  as  recommended  by  Ruff  in  1926  leads  to  log^Q  760  mm  = 

2.  88  at  2,  590*  C  as  he  had  reported  at  that  time. 

If  we  plot  the  three  calculated  vapor  pressure  figures  as  given  by  Ruff  in 
1935  and  draw  the  dashed  line  shown  in  the  graph  (figure  12),  we  obtain  an 
intercept  of  8.  15  at  1  /T  =  0. 

Since 

Pmm  * 


and 


lo«lOP«nm  ■  loI  10  P.on  +  lo*  10760  . 


the  intercepts  for  the  different  pressure  units  are  therefore 


{toflOPmm]  1/r  _  o  “  U’°5 
he*-]  l/T  -  O  "  8-15 


A  value  of  8.  15  for  the  intercept  at  l/T  =  0  would  appear  reasonable  if  the 
ordinate  were  given  as  log10  Patm  instead  of  the  log1Q  Pmm  actually  used. 

It  thus  appears  that  Ruff  may  have  mistakenly  used  8.  15  as  an  intercept  when 
log10  Pmm  was  plotted  as  the  ordinate.  Actually,  to  be  correct,  it  is  neces¬ 
sary  to  plot  logjQ  Patm  in  conjunction  with  the  8.  15  intercept. 

In  view  of  the  uncertainties  in  the  analysis  just  given  above,  it  is  desirable 
to  present  some  independent  data  supporting  Ruff's  1926  extrapolation  to  an 
intercept  of  11.  1  at  l/T*  0.  We  have  already  shown  in  Section  II,  E  that 
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Pt  (atm)  -  3/2 e 


18.41  (1  -  ii.60) 

T 


or 


Pt  (mm)  -  3/2*760  e 


18.41  (1  - 


The  second  equation  can  be  rearranged  to  the  form 
lo8lOPt(nun)  "  11*1  ~  25,200/T 


in  which  it  is  obvious  that  the  intercept  is  11.  1. 

We  therefore  conclude  that  the  boiling  point  of  2,  950*  C  given  by  Ruff  in  1935 
is  too  high.  On  the  other  hand  the  value  of  2,  590*  C  is  only  a  lower  limit. 

This  is  particularly  true  since  his  measurements  were  performed  under  re¬ 
ducing  conditions. 

Information  on  the  very  high  boiling  point  of  3,  500#  C  given  by  Mott51  in  1918 
indicates  that  this  measurement  was  a  very  rough  one.  SiOz  was  decomposed 

in  an  arc,  and  the  temperature  measurements  were  necessarily  quite  un¬ 
certain. 

Summarizing  then,  the  recorded  boiling  points  for  SiO^,  in  the  literature, 
are  spread  over  a  considerable  range.  The  lower  values  of  2,  230*  C  and 
2,  590*  C  must  be  considered  as  lower  limits.  The  value  of  2,  950*  C  is 
apparently  a  mistake  and  the  value  of  3,  500*  C  is  only  an  approximation. 

Thus,  our  calculated  value  of  2,  800*  C  appears  to  be  reasonable  and  is  prob¬ 
ably  more  reliable  than  the  values  in  the  literature* 

Additional  vapor  pressure  data  for  the  decomposition  of  Si02,  under  so-called 
neutral  conditions,  have  been  given  by  Brewer  and  Mastick2,  Inuzuka  and 
Ageha  ,  and  Porter,  Chupka,  and  Inghram  3  .  The  data  of  Brewer  and 
Mastick  have  been  shown  to  agree  with  the  calculations  based  on  the  first 
reaction  SiO^siO^  +  *02(g)  .  inasmuch  as  we  have  found  that  con¬ 

sideration  of  all  the  other  possible  reactions  only  changed  the  predicted  boil¬ 
ing  point  from  3,  085*  K  to  3,  070*  K,  it  can  be  seen  that  the  first  reaction  is 
of  primary  importance.  Our  calculated  vapor  pressure  data  at  low  temperatures 
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are  compared  with  the  data  of  Brewer  and  Mae  tick2  in  figure  5.  The  agreement 
ie  good  and,  if  anything.  Brewer  and  Mastick's2  extrapolated  data  might  indi¬ 
cate  an  even  higher  boiling  point  than  the  one  we  have  calculated. 

The  data  of  Porter  et  al3  are  seen  to  be  too  high  relative  to  our  vapor  pres¬ 
sures  infigure  5.  However,  they  were  apparently  not  operating  under  strictly 
neutral  conditions.  Thus,  for  1,  900*  K  they  gave  pc.0  -  3.4  *  10  —  5  and  p0  - 
8.2  z  10  ™  ® atm.  It  follows  thatP^Q  Pq^  -9.7  z  10  -  9  atm  3/z\  2 

and  that 


AFPofter,  1,900  "  ♦  69,595  calorie/mole 


As  seen  in  Table  VH-A,  this  free  energy  value  agrees  closely  with  our 
thermochemical  value  of  71,000  calories /mole. 

The  data  of  Inuzuka  and  Ageha3®  were  discussed  earlier  and  it  was  shown 
that  their  reported  high  vapor  pressure  data  could  be  accounted  for  by  a  re¬ 
ducing  reaction  involving  the  molybdenum  wire.  In  fact,  the  excellent  agree¬ 
ment  obtained  with  the  thermochemical  data  for  SiO^  suggests  that  these 

data  are  fairly  reliable. 

B.  OXIDIZING  CONDITIONS 

Just  as  in  the  case  of  neutral  conditions  we  consider  the  same  two  reactions 
for  the  case  of  oxidizing  conditions  in  air.  That  is,  we  neglect  other  possible 
reactions.  It  is  possible  that  excess  oxygen  or  nitrogen  might  lead  to  new 
species  which  are  not  considered,  here,  but  until  experimental  data  clarify 
this  point  we  make  the  assumption  that  other  reactions  are  negligible. 


si02(l7“**si0(g)  +  *°2<g)  ,  (17) 

X°XtT~°(t>  ■  <18> 


It  is  assumed  that  SiC^  decomposes  in  the  presence  of  air  at  various  total 
pressures,  p*.  For  the  present  case,  nitrogen  is  assumed  to  be  inert.  Then 
we  have  the  following  relations: 
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pc  •  pN2  ♦  PQ2  ♦  PSiO  ♦  Po  • 


(19) 


Pn2"4<Po2  ♦  HPq  “  ^PsiO>  • 


(20) 


K1  "  PSiO  P02 


(21) 


-  Pc/POj  * 


(22) 


We  can  solve  the  above  set  of  four  equations  for  the  four  unknowns  pN  #  pQo. 
and  Pgjo  in  terms  of  Pt,  K|,  and  *2  as  known  parameters.  Combining  * 

equations  (19)  and  (20)  we  get 


pt  "  5Pq2  +  3Pq  -  Psio  •  (23) 


But  from  equation  (22),  pQ  •  K2Po^  #  and  equation  (23)  becomes 
pt  -  ’  POj  +  3  KjPo'^  -  PsiO 

or 

PsiO  "  5  Po2  +  BKjPo^  -  Pt  • 

Also  from  equation  (21)  we  have 


(24) 


(25) 


(26) 


Combining  this  with  equation  (2  5),  we  finally  get 
5  ♦  3K2p02  -  PtPo«  -  Kx  -  O 


(27) 
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After  this  cubic  equation  has  been  solved  for  pQ  ,  the  other  unknowns  can 

2 

be  evaluated  from  equations  (26),  (22)  and  (20). 

Solutions  for  this  set  of  equations  were  obtained  for  values  of  pt  =  0.  1,  1,  10, 
and  50  atm.  Calculations  were  made  in  the  range  of  2,  000*  to  3,  750*  K.  The 
resulting  data  have  been  plotted  in  figures  13,  14,  15,  and  16  which  show  the 
partial  pressures  of  SiO,  O2,  Si02»  and  N2  as  functions  of  1/T. 

It  is  seen  that  at  temperatures  well  below  the  boiling  point,  the  Si0^gj  forma¬ 
tion  is  markedly  repressed.  Thus  from  figure  13atpt  =  0.  1  atm  and  T  =  2,  000*  K, 
the  SiO  partial  pressure  is  essentially  equal  to  the  Si02  partial  pressure.  How¬ 
ever,  in  this  case  a  boiling  point  is  not  reached  until  about  2,  700*  K.  At  all 
higher  temperatures,  SiO(g)  formation  becomes  more  important  than  Si02^gj 

formation.  For  reference  purposes  in  figures  13  through  16.  the  vaporization 
curve  for  SiO  under  neutral  conditions  is  also  given.  The  curves  given  in 
each  of  the  figures  extend  up  to  the  boiling  point.  At  higher  temperatures, 
total  pressures  are  obtained  in  excess  of  the  postulated  pressure.  At  the 
boiling  point  and  higher  temperatures  the  N2  pressure  is  zero.  Also  at  the 
higher  temperatures,  the  SiO,  O2  and  O  partial  pressures  are  the  same  as 
for  tne  neutral  conditions  already  discussed.  However,  these  data  are  not 
repeated  in  figures  13  through  16.  The  rapid  variation  in  partial  pressure  of 
N2  has  not  been  indicated  with  much  precision  but  an  approximate  indication 
is  given  in  the  figures. 

The  data  presented  in  figures  13  through  16  show  that  under  some  conditions 
the  oxygen  present  in  air  is  sufficient  to  repress  SiO  formation  significantly. 
Although  the  effect  decreases  near  the  boiling  point  it  becomes  very  pro¬ 
nounced  at  lower  temperatures. 

Consequently,  in  any  proposed  application  of  fused  silica,  the  conditions  of 
the  high -temperature  exposure  must  be  known  before  the  equilibrium  gas 
species  can  be  estimated. 
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APPENDIX  A 


SUMMARY  OF  PHYSICAL  AND  CHEMICAL  PROPERTIES  OF  SILICA 


Because  of  the  great  current  interest  in  silica  as  a  possible  consumable  re- 
entry  vehicle  material,  there  is  a  need  for  a  convenient  summary  of  the  best 
available  data  on  its  properties.  The  following  summary  presents  these  data 
in  a  concise  systematic  form.  It  contains  not  only  the  material  from  the 
literature,  reviewed  in  the  main  body  of  the  report  plus  the  results  of  calcula  - 
tions  with  the  literature  data,  but  it  also  includes  other  data  on  mechanical 
and  heat  transfer  properties,  thermal  shock  resistance,  crystal  structure, 
electrical  and  magnetic  properties,  availability,  and  cost. 

A.  DECOMPOSITION  REAC  TIONS  AND  HEAT  OF  VAPORIZATION 

The  two  processes  below  have  been  considered  but  under  equilibrium  condi¬ 
tions  (Eq.  1)  decomposition  to  form  SiO  and  O £  has  been  shown  to  be  the 
most  important  reaction. 

si°2(irsi°(g)  +  ^°2(g)  AHv«poriz«tion,  3,000°K  "  171,000  cslorie/mole  (1) 

-  5,130  BTU/lb 

Si02(ir^Si02(g)  AHT«porlmtion,  3,000°K  -  127,000  cslorie/mole  (2) 

-  3,810  BTU/lb. 


B.  VAPOR  PRESSURE 

The  vapor  pressure  of  SiO^  is  not  a  unique  function  of  temperature  but  is  also 
very  sensitive  to  the  presence  of  oxidizing  or  reducing  agents  in  the  atmos¬ 
phere.  The  total  vapor  pressure  (pt)  versus  temperature  (T)  for  silica  under 
neutral  conditions,  where  the  decomposition  is  essentially 

si02arsl0(g)  *  h°2(,>  . 

is 

,  18.41  fl-  3160  |  ,  18.41(1  -5680/TR) 

Pt  -  3/2e  \  T  /  -  3/2e  R  ’ 
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where 

P*  "  PsiO  ♦  Po2  • 
pSiO  -  2/3  Pt  , 

pG2  -  1/3  Pt  . 

pt  •  PSiO  »  ““d  Po2  ■**  in  atmo^herec, 
T  -  °Kelvin  (°K) 

Tg  -  °Rankine  (°R) 


This  formula  is  estimated  to  be  accurate  to  ±  50  percent  near  3,  000*  K  and 
applies  to  the  temperature  range  above  2,  000*  K.  It  predicts  a  total  vapor 
pressure  of  2.  17  x  10"  atm.  (1.  65  mm  Hg)  at  2,  060*  C  (2,  330*  K).  This 
is  rmu:h  lower  than  the  value  of  1 1  mm  reported  at  this  temperature  by  Ruff 

et  al.  However,  their  measurement  was  per  formed  only  once  and  they  stated 
that  reducing  conditions  made  their  vapor  pressure  determination  give  high 
values.  For  approximate  estimates  of  PsiQ  under  neutral  conditions  see  Fig.  4. 

Small  contributions  to  the  vapor  pressure  are  made  by  other  equilibrium  prod¬ 
ucts.  (See  Fig.  11).  The  following  is  a  typical  expected  neutral  equilibrium  mix¬ 
ture  at  a  temperature  of  3.  000*  K  (2,  727*  C,  4,  940*  F,  5,  400*  R): 


Species 

Partial  Pressure  (atm 

SiO 

0.  43 

°2 

0.  18 

O 

0.  05 

Si02 

0.  04 

Si202 

5  x  10"5 

Si 

1.  3  x  10”^ 

Pt  =  0.  70  atm 

For  the  effects  of  an  oxidizing  atmosphere  under  different  total  pressures 
see  Figures  13-16. 
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G.  MOLECULAR  WEIGHT 


The  molecular  weight  is  60.  06. 

D  MEETING  POINT 

t  used  silica,  being  a  supercooled  liquid,  has  no  recognized  melting  point. 

:  .softms  (viscosity  S'  107-  6  poises)  at  1.  500*  C52.  The  various  crystalline 
torms  of  SiO^  do  melt,  however. 

F o rm  Melting  point 


°  F 

*  C 

c  ri  stobalite 

3,  140 

1. 728 

tridymite 

3,  056 

1, 680 

qua  rtz 

2,  930 

1,610 

E.  BOILING  POINT 

As  indicated  in  (A)  and  (B)  of  this  appendix,  vaporization  of  SiOz  yields  pre 
dominantly  a  mixture  of  SiO  and  Oz.  Lesser  amounts  of  SiOz  gas  and  O 
atoms  are  formed.  The  calculated  normal  boiling  point  (vapor  pressure  = 

1  atm)  is  5,072*  F  (2,  800°  C.  3,  070*  K). 

F.  VISCOSITY 

A  viscosity  relation  which  is  valid  for  the  temperature  range  of  2,  197*  K  to 
2,  333*  K  (1 , 924*  C  to  2,  060*  C)  is 

io«l0^  ■  28,200/T(OK)  _  7.73  53‘  54 


G.  SPECIFIC  HEAT 

For  vitreous  or  fused  silica  the  molar  heat  capacity  relation  given  by  Kelley1-* 
covers  the  range  of  298*  K  to  2,  000*  K.  Its  accuracy  is  ±2  percent. 

<~p  -  13-38  +  3-68  x  10  ~3T  -  3-45  *  105T  “  2 
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where 


Cp  ■  calorie/mole  —  degree  K 
T  -  °KeUin 


H.  SENSIBLE  HEAT  PLUS  HEAT  OF  VAPORIZATION 

The  sensible  heat  required  to  bring  fused  silica  from  298*  K  to  3,  000*  K  is 
found  to  be  1,  524  Btu/lb  (or  50,670  calorie /mole).  Going  from  298*  K  to 
3,  070*  K  requires  1,570  Btu/lb  (or  52,  190  calorie/mole).  Assuming  that 
the  heat  of  vaporization  at  3,  070*  K  does  not  differ  appreciably  from  that  at 
3,  000*  K,  one  can  obtain  the  following  total  heat  absorptions  for  the  principal 
decomposition  equations  (A-l)  and  (A- 2)  of  this  appendix. 


Temperature 


Process 


Btu/lb 


Calorie /mole 


3, 000*  K 
3.  000 
3.  070 
3.  070 


Reaction  ( 1 ) 

6,  650 

222, 183 

U) 

5,  336 

177, 877 

(1) 

6.  700 

U> 

5,  386 

915  mole  8  SiO  by  (1) 

6,  790 

0.11  moles  O  by  1/2  —  »  O 

0.  085  moles  SiO^  by  (2) 


These  data  show  that  the  total  heat  of  decomposition  is  essentially  controlled 
by  reaction  (1 ). 


That  is  even  when  the  equilibrium  mixture  at  3,  070*  K  is  considered  the  extra 
energy  required  to  dissociate  is  partially  offset  by  the  reduced  heat  r e - 
quired  to  evaporate  gaseous  SiO^.  It  can  also  be  seen  that  the  uncertainty 
in  the  heat  of  reaction  (1)  at  3,000*  K  was  of  the  order  of  5,  000  calorie/mole 
or  150  Btu/lb.  Thus,  for  most  purposes,  it  is  sufficient  to  use  only  the  data 
from  reaction  (1). 

I.  DENSITY 

The  density  of  pure  fused  silica  is  1  37.  3  lb/ft3  or  2.  20  gm/cm3  .  52 

J.  THERMAL  CONDUCTIVITY 


At  O*  C,  fused  silica  has  a  thermal  conductivity  of  0.  795  Btu/hr  ft*  F  or 
10“  calorie  -  cm/sec  -  c -  *  C. 


33  x 
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K.  COEFFICIENT  OF  THERMAL  EXPANSION 


The  coefficient  of  thermal  expansion  for  clear,  fused  silica  from  20*  -  320*  C 
is  reported  to  be  3.  05  x  10“^/*  F  or  0.  55  x  10”^/*  C  by  G.  E,  and  from  ^ 
0°  -  300*  C  it  is  reported  to  be  3.  1  x  10"^/*  F  or  0.  56  x  10*6/*  C  by  Corning. 

L.  YOUNG'S  MODULUS 


Clear,  fused  silica. 
Modulus  of  10.  1  x  10 


at  room-temperature  is  reported  to  have  a  Young 
^  lb/in^  or  7,  138  kg/mm^  by  G.  E.  ^ 


's 


M.  TENSILE  STRENGTH 

In  1/4  inch  rods,  fused  silica  has  a  strength  of  7,  100  lb/in^  or  5  kg/mm^; 
whereas  in  fibres,  it  has  a  strength  of  43,  500  to  116,  000  lb/in^  or  30.  6  - 
8i .  6  kg/mm^  . 

N.  THERMAL  SHOCK  RESISTANCE 

Fused  silica  has  excellent  thermal  shock  resistance.^ 


Other  forms  of  silica  have  poor  resistance  at  certain  temperatures  because 
of  volume  changes  which  occur  at  inversion  points. 

A  6  x  6  x  1  /4-inch  annealed  plate  of  fused  silica  can  withstand  a  1, 000*  C"^ 
differential  from  being  plunged  into  cold  water  after  oven  heating. 

O.  HARDNESS 

55 

The  hardness  of  fused  silica  is  4.  9  on  the  Mohs  Scale  whereas  cristobalite, 

4  s 

quartz,  tridymite  have  a  value  of  7.  0  on  the  Mohs  Scale. 

P.  CRYSTALLINE  STRUCTURE 

Pure  SiO^  is  a  polymorphic  material.  The  vitreous  or  glassy  form  of  SiO^  is 
a  metastable,  undercooled  liquid  at  all  temperatures  below  1,728*  C  (2,000*  K). 
The  rate  of  transformation  of  vitreous  SiO^  into  the  other  forms  may  become 
measurable  above  1,000*  C.  Impurities  also  affect  the  rate  of  transformation. 
Listed  below  are  the  stability  ranges  for  the  various  forms  of  SiO^,: 
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F  orm 


Temperature  range  of  stability 


Silica  glass 

T  >  1, 728*  C 

&  -  Cristobalite 

1.728*  -  1.470*C12 

/9  -  Tridymite 

1, 470*  -  870*  C 1 1 

/9  -  Quartz 

870*  -  573*  C 

a  -  Quartz 

573*  - 

In  fused  silica  glass,  the  Warren  model  is  generally  accepted**^  to  represent 
the  short  range  order.  Being  a  liquid-like  material  there  is  no  long-range 
order  greater  than  about  4  or  5  interatomic  distances.  According  to  the 
War.  en  view  each  silicon  atom  is  surrounded  tetrahedr ally  by  four  oxygen 
atoms.  Each  oxygen  atom  is  shared  by  two  silicon  atoms.  The  SiO  distance 
is  about  1.6  A,  the  0-0  distance  is  2.65  A,  and  the  Si-Si  distance  is  3.2  A. 

Q.  ELECTRICAL  RESISTIVITY  1  1 


The  electrical  resib ..vities  of  various  forms  of  silica  are: 


Form 


Temperature 


Vitreous  20*  C 

600*  C 
1, 300*  C 

Crystalline  20*  C 

(quartz)  1 , 300*  C 

20*  C 
1, 300*  C 

R.  DIELECTRIC  CONSTANT1  1 
The  dielectric  constants  of  silica  are: 


Resis  tivity 
(ohm-cm) 


X 

o 

*— 

<J1 

x  106 

x  105 

x  1014 

II  to  prin¬ 

x  10-* 

cipal  axis 

x  1016 

j.to  prin¬ 

x  104 

cipal  axis 
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Form 


Tempe  ratur  e 


Dielectric  Constant 


Vitreous  20*  C  3.  75 

Crystalline  (quartz)  20*  C  4.6  ||  to  direction  of 

field 

4.  5  J_to  direction  of 
field 

S.  MAGNETIC  SUSCEPTIBILITY  1 1 
The  magnetic  susceptibilities  of  silica  are: 

Sus  ceptibility 

-0.45  x  10“^  magnetic  units  /gm 
-0.  45  x  10"^  magnetic  units /gm 


Form 
V  itreous 

Crystalline  (quartz) 


T.  PURITY 


Fused  silica  can  be  obtained  as  one  of  the  purest  chemicals  known.  For 
example.  Coming's  commercial  grade  is  99.9  percent  pure.  It  is  possible 
to  obtain  a  grade  of  99.  99  percent  purity  for  use  in  delay  lines,  and  a  grade 
of  99.999  percent  purity  is  available  for  use  in  spectroscopy.  The  cost 
climbs  rapidly  with  purity. 

Some  analyses  for  the  Corning  grades  are  presented  below.  It  should  be 
emphasized  that  the  impurities  given  need  not  be  unavoidable  since  precau¬ 
tions  can  be  taken  to  remove  certain  impurities  if  necessary.  The  type  7943 
grade  (99.  99  percent  purity)  is  a  desirable  optical  grade  for  infrared  work 
because  of  the  absence  of  the  OH  group.  As  of  March  1958  it  was  not  avail¬ 
able  in  sizes  over  1 /8-inch  thick.  Type  7940  can  be  obtained  in  flat  sheets 
to  20  inches  in  diameter  and  up  to  3  inches  thick.  ^8 


Impurity 

Type  7940 
(pe  rcent) 

OH 

0.  1 

Li 

<  0.  0001 

Na 

0.  00003 

Type  7943 
(percent) 

0.  00 


0.  001 
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Impurity 

Type  7940 
(percent) 

Type  7943 
(percent) 

K 

0.  0001 

0.  0003 

Fe 

0.  002 

0.  001 

Cu 

0.  00001 

0.  00001 

Ti 

0.  0005 

0.  0005 

Bi 

0. 0000047 

P 

0.  000010 

As 

0. 00000057 

Sb 

0. 00000051 

Zr,  V,  Pb,  Sn,  Ge, 

Mo,  Cs 

not  detected,  and 
less  than  ~  1 0  ppm  = 
0. 001  percent 

A12°3 

0.001  59 

typical  analysis  of  some  G.  E.  fused  silica  is  presented  below.  The  source 
a  Brazilian  quartz  sand.  The  data  were  given  to  D.  Walker  by  K.  Field 
General  Electric.  °° 

Impurity 

Percent 

Fe2°3 

0.  0017 

A12°3 

0.  014 

Ti02 

nil 

Z*°2 

nil 

K 

0.  0003 

CsO 

<  0.  0003 

MgO 

<  0. 0007 

C 

B 

none 

<  0. 0001 
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U.  AVAILABILITY 


Fused  or  vitreous  silica  is  produced  by  at  least  three  manufacturers  in  this 

country : 


a.  Corning  Glass  Works,  Corning,  New  York. 

b.  General  Electric  Company,  Lamp  Glass  Department,  Nela  Park, 
Cleveland  1 2,  Ohio. 

c.  Amersil  Division  of  Engelhard  Industries. 

v.  Cost 

The  cost  of  fused  silica  is  very  dependent  on  the  purity,  finished  form,  etc. 
Some  typical  prices  for  G.  E.  products  are  shown  below: 


Form 

Clear,  Fused  Quartz 
(Cost/lb.  ) 

T  ranslucent 
(Cost/lb.  ) 

rod,  1  / 1 6-inch  diameter 

$50. 00 

-  - 

rod,  3/16-inch  diameter 

30.  00 

$12. 00 

ingot,  to  18  1/2  inches 
diamete  r 

15.  00 
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APPENDIX  B 
ILLUSTRATIONS 
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FREE  ENERGY  CHANGE,  AF,  K  CAL 


Figure  1  FREE  ENERGY  CHANGES  OF  REACTIONS  IN  Si02  DECOMPOSITION 

VERSUS  TEMPERATURE 
4508 
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EQUILIBRIUM  CONSTANT 


TEMPERATURE, • K 


K)-« 


•  w  .«• 

tT*k)  x  10 


Figure  2  EQUILIBRIUM  CONSTANTS  OF  Sip2  DECOMPOSITION  REA( 

VERSUS  « 

4509  T 
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y 


TEMPERATURE  -  *K 


Figure  4  PARTIAL  PRESSURE  OF  SiC>(  j  OVER  SiO,  VE 
TEMPERATURE  NEGLECTING  Oz  DISSC 

4511 


VERSUS  RECIPROCAL  OF 
'SOCIATION 


83 


PRESSURE  ,  (ATM  ) 
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5  PARTIAL  PRESSURE  OF  SiO<B)  OVER  CRISTOBAUTE  VERSUS 
TEMPER  AT  URE  FROM  1 , 800  *  TO  2, 000  *  K 


I 

t 


I 
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Figure  7  DISSOCIATION  PRESSURE  RATIOS  OF  OXYGEN  VERSUS  TEMPERATURE 

NEAR  3, 000* K 
4514 


I 

1 

! 
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Figure  9  PARTIAL  PRESSURES  OF  IMPORTANT  SPECIES  WHEN  SiO 
DECOMPOSED  BY  CARBON  OR  SILICON 
4516 


PRESSURE (ATM) 


oo  O  O  O  O  O  o 

<r>  o  o  q  o  o  o  § 

S5  v  ir,  <\J  _ 
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-  TEMPERATURE,#K 

8  §  8  8  8  8  8 

<7>  00  h-  %D  *D  «  K) 

CNJ  OJ  CNJ  OJ  rg  <\j  fsj 

i._ir  i  i  — f-  f-= 


1  •  p total 
n  •  ps.o 
in  ■  po2 
nr  ■  po 

T  -PS.02 

^■ps,2o2 

9  •action*  Con*d«r«d 
'  s  1  °2(jt)  ““  Si0(9)  +  I/2  02{q) 

2  1/  2  02  (g  j  ~ 0(g) 

3  Si02(1)  —  Si02(Q  )  T  » 

4  2SlO(g)-»Si202(g) 


T  »  2000*  K 


10'*  :  O  EXPERIMENTAL  POINTS  OF 
PORTER  ETAll3,F0R  PSl0z(g) 

x  EXPERIMENTAL  POINTS  OF 
BREWER  AND  MASTICK121 
•O’,:  F0R  PSI0,«, 


I0'*L 
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Figure  11  PARTIAL  PRESSURES  OF  PRODUCTS  WHEN  SiO^  DECOMPOSES  IN 
AIR  VERSUS  ^  UNDER  NEUTRAL  CONDITIONS  A^ID  INCLUDING 

OXYGEN  DISSOCIATION 
45  ] 

_  Q 


f*®**K»  /» 30* 

\2t9o*c/  (toao*c 


PRESSURE  (ATM) 


TEMPERATURE,  •  K 


Figure  13  PARTIAL  PRESSURES  OF  PRODUCTS  WHEN  SiO^  DECOMPOSES  IN 

AIR  AT  0.  1  ATM  VERSUS  2. 
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Figure  15  PARTIAL  PRESSURES  OF  PRODUCTS  WHEN  SiOz  DECOMPOSES  IN 

AIR  AT  10  ATM  VERSUS  i 
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